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In  this  project  a  theoretical  model  was  formulated  to  predict  the 

optical  properties  of  non-homogeneous  materials  (cermets).  This  task  has  now 
been  fully  completed.  Our  calculations  and  the  model  have  corroborated  our 
original  prediction  that  correlation  effects  due  to  the  formation  process  of 
cermets,  quantum  size  corrections,  and  renormalized  corrections,  are  essential 
ingredients  to  give  a  realistic  theory  for  the  optical  properties  of  cermets. 
This  is  particularly  true  in  the  high  concentration  regime  of  the  metallic 
minority  constituent  (more  than  10%  volume  fraction  of  the  metallic 

constituent).  The  results  of  this  work  have  been  published  in  six  papers. (See 
Appendix  I). 

The  principal  findings  which  are  presented  and  elaborated  in  the  published 
articles  can  be  summarized  as  follows: 

1.  For  cermets  with  island  radiuses  under  100A,  the  dielectric  constant  of 

the  metallic  constituent  must  be  calculated  using  a  quantum  formalism. 

For  radiuses  that  are  between  100  and  50A  the  formalism  developed  by  D. 

M.  Wood  and  N.  W.  Ashcroft  and  the  one  developed  by  A.  Kawabata  and  R. 
Kubo  give  similar  results  except  for  a  geometrical  factor,  while  for 
islands  of  radius  under  50A,  the  Wood  and  Ashcroft  formalism  must  be 
used.  The  structure  found  in  the  Wood  and  Ashcroft  model  must  be 
smoothed  by  averaging  the  dielectric  constant  with  a  log-normal 
distribution  of  radiuses  before*  it  can  be  used  in  calculating  the 
dielectric  properties  of  cermets. 

2.  Using  the  Random  Continuum  Approach  and  doing  a  multiple  scattering 
calculation  in  the  First  Smoothing  Approximation,  a  dielectric  constant 
for  cermets  was  obtained.  The  resulting  dielectric  constant  improved 
the  zeroth-order  non-correlated  model  of  Bruggeman  by  reducing  the 
intensity  of  the  resonance  of  the  complex  index  refraction,  thus, 
achieving  a  better  agreement  with  experimental  results. 


3.  A  better  effective  dielectric  constant  was  obtained  when  the  First 
Smoothing  Approximation  was  improved  by  developing  a  renormalized 
procedure  that  used  a  T-Matrix  formalism.  With  the  renormalized 

approach,  the  resonance  peak  of  the  refractive  index  was  shifted  to 
the  correct  position  in  the  spectrum,  thus  obtaining  excellent 
agreement  with  the  experimental  results.  This  calculation  was 

performed  taking  MG  as  the  zeroth-order  approximation  dielectric 

constant  since  it  is  based  on  a  model  that  sees  cermets  as  distinct 

island  embedded  in  a  dielectric  medium. 
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QUANTUM  SIZE  EFFECTS 
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Physics  Department,  University  of  Puerto  Rico,  Rio  Piedras, 

Puerto  Rico;  tPresent  address:  Cornell  University,  Ithaca,  NY 

Abstract  The  dielectric  constant  of  metallic  grains  in  real 
cermets  is  obtained  using  the  formalism  developed  by  Wood  and 
Ashcroft,  and  the  resonance  structure  derived  from  this 
formalism  is  smoothed  out  by  using  a  log-normal  distribution 
to  average  the  dielectric  constant  over  grain  sizes.  The 
results  from  the  quantum  theory  for  the  metallic  inclusions 
in  AU/AI2O3  cermet  are  compared  with  Drude's  results  against 
experimental  data  using  mean  field  theory  to  obtain  the 
effective  medium  dielectric  constant  of  the  cermet. 


INTRODUCTION 

Cermets  are  special  types  of  composite  materials  made  of  two 
constituents,  one  metallic  and  the  other  dielectric,  possessing  a 
granular  structure  where  the  typical  diameter  of  the  grains  is  of 
the  order  of  a  100  A.  Given  the  small  size  of  the  microstructure, 
it  is  always  possible  to  describe  its  infrared  and  visible  optical 
properties  in  terms  of  an  effective  medium  dielectric  constant. 

Mean  field  theories  that  express  the  effective  dielectric 
constant  of  the  cermet  in  terms  of  the  known  parameters  of  the 

constituents  have  been  developed.  The  best  known  effective  medium 

.  1  2 
theories  are  the  Maxwell-Garnett  and  the  Bruggeman  theories. 

Most  of  the  mean  field  theories  perform  calculations  using  for  the 

metal  constituent  a  bulk  dielectric  constant  that  is  corrected  by 

an  effective  mean  free  path  due  to  surface  scattering  of  the 

electrons  at  the  grain's  boundaries. 

Our  objective  is  to  develop  a  more  realistic  dielectric 

constant  for  the  metallic  grains*-  that  take  into  account  quantum 

size  effects  (QSE).  The  quantum  corrected  dielectric  constant  of 

the  metal  will  be  susbtituted  into  the  Bruggeman  theory  and  the 
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results  compared  with  experimental  data  for  real  cermets. 

This  work  is  presented  in  two  parts.  First  part  describes 
the  theoretical  model  used  to  consider  quantum  effects.  Second 
part  (following  paper)  describes  the  procedure  to  smooth  out  the 
discrete  structure  appearing  in  the  QSE  calculations  and  compares 
the  effective  medium  dielectric  constant  of  the  Au/A1203  cermet 
calculated  using  QSE  and  Drude's  model  (D)  against  experimental 
data.  Final  conclusions  appear  in  the  second  article. 

Quantum  Size  Effects 

Due  to  the  sizes  of  the  metallic  grains  the  bulk  dielectric 
constant  of  the  metal  must  be  corrected  for  size  effects  before 
being  used  in  the  mean  field  calculations.  Traditionally  the  D 
constant  together  with  the  interband  terms  is  used  and  the  mean 
free  path  is  corrected  for  grain  boundary  collisions.  Unfor¬ 
tunately  this  correction  is  not  justified  for  smaller  particles 
since  it  does  not  fully  take  into  account  size  effects. 

For  islands  of  sizes  smaller  than  100  A  a  quantum  treatment 

3 

is  needed.  We  will  introduce  the  Wood  and  Ashcroft  formulation 
in  calculating  quantum  size  effects  for  small  metallic  particles. 
We  summarize  the  salient  points  of  their  theory.  The  induced 
current  produced  by  the  applied  field  is  calculated  expanding  the 
density  matrix  in  terms  of  the  perturbing  field.  The  Liouville 
equation  is  then  solved  using  the  Linear  Response  Theory  in  the 
Random  Phase  Approximation,  and  then  simplified  within  the 
Diagonal  Response  Approximation  that  is  equivalent  to  assume  that 
the  electron  system  is  homogeneuos,  and  only  responds  to  the  wave 
vector  of  the  applied  field. 

In  the  region  of  interest  for  the  optical  properties  of 
cermets  (0 . 2pm<\<20pm)  the  long  wave  length  limit  applies  and  it 
is  possible  to  expand  the  expression  for  the  dielectric  constant 
obtained  from  the  induced  current  in  powers  of  the  ratio  between 
the  frequency  of  the  applied  field  and  the  plasma  frequency.  This 
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yields  a  dielectric  constant  expression  in  terms  of  the  matrix 

elements  of  the  momentum  of  the  electron  along  the  applied  field. 

The  grains  in  cermets  usually  have  completely  undefined 

shapes  but  since  the  wave  lengths  of  interest  to  this  work  are 

considerably  larger  than  the  grain  sizes,  we  can  again  follow  the 

approach  of  Wood  and  Ashcroft  and  assume  the  grains  to  be  cubic  in 

form.  The  field  is  then  taken  to  lie  in  the  direction  normal  to 

one  of  the  box  faces.  Then  all  that  is  needed  is  the  spectrum  of 

an  electron  in  a  one-dimensional  quantum  box.  To  take  into 

account  the  effect  of  a  finite  electron  relaxation  times,  the 

adiabatic  switching-on  term  is  replaced  by  an  effective  relaxation 

4 

time  T.  Within  the  Mermin  scheme  the  dielectric  function  for  the 
metallic  particles  can  then  be  written  as, 

m 

Ree(x)  =  1  +  44  a  2  m2  (m2  -  m2)  x 
n*  a0  m=l 

“  m,2(A2-(x2+f~2)  ]  [  l-(-l)m^m'  ] 

m'  =  l  0 

m 

Im£ (x)  =  44  a  T  2  m2  (m2  -  m2)  x 

— X  -  “ 

7T’  a0  x  m=l 


“  m'2[A2+x2+r2] (l-C-l)™*™'] 

m'=l  0  (1) 

,  where 

0  =  A3  {(A2  -  x2  +  T2)2  4x2  r2};  A  =  m'2  -  m2 
x  =  2huj  ma2/h2n2;  T  =  hmc2/T£f;  mc  =  INT(Kfa/n) 

£ j  is  the  Fermi  energy,  the  Fermi  wave  number,  a  the  size  of 
the  grains,  and  ao  the  Bohr  radius. 

Figure  1  compares  the  dielectric  constant  calculated  using 
the  QSE  as  expressed  in  Eq.  ()),  with  that  calculated  using  D 
model.  In  the  calculation  of  the  D  model  the  size  of  the  silver 


::«i  |5M>| 


M.  GOMEZ.  G.  RODRIGUEZ  AND  L  FONSECA 


FIGURE  1  Dielectric  constant  of  Ag  metallic  grains 
calculated  using  D  model  and  QSE  for  a  30  A  grain. 

particles  are  taken  to  be  30  A  and  no  interband  contributions  have 
been  included.  D  model  was  corrected  for  size  effect  using  an 
effective  relaxation  time  that  corrects  the  bulk  result  with  grain 
boundary  scattering.  The  multiple  peaks  appearing  in  the  QSE 
calculations  are  due  to  the  discrete  energy  spectrum  of  an 
electron  in  a  box.  In  the  next  part  of  this  work  (following 
article)  we  will  discuss  the  application  of  QSE  to  real  cermets. 
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Abstract  See  Part  I  (preceding  paper). 

Application  of  the  QSE  to  a  Real  Cermet 


We  will  apply  the  Quantum  Size  Effects  (QSE)  result  described  in 
Part  I  to  the  calculation  of  the  optical  properties  of  real 
cermets  and  compare  the  results  with  those  obtained  using  Drude's 
model.  In  these  calculations  interband  contributions  to  the 
dielectric  constant  of  the  metal  grains  will  be  fully  taken  into 
account.  Experimentally  the  resonances  observed  in  the  dielectric 
constant  of  an  isolated  metal  grain  due  to  QSE  are  not  observed  in 
the  effective  dielectric  constant  of  a  real  cermet.  The  disap¬ 
pearance  of  these  peaks  is  due  to  the  fact  that  the  metallic 
grains  in  a  real  cermet  do  not  have  a  unique  size,  but  instead 
have  a  distribution  of  sizes.  For  this  reason  the  dielectric 
constant  obtained  from  the  QSE  calculations  for  one  isolated 
particle  was  averaged  over  a  distribution  that  realistically 
represented  variation  in  size  in  actual  cermets. 

The  histogram  of  the  size  distribution  is  skewed  and 
actually,  a  log-normal  distribution  is  usually  used  to  fit  the 
experimental  results*.  Averaging  the  dielectric  constant  of  one 
metallic  grain  with  a  log-normal  distribution  function  an 
effective  dielectric  constant  for  the  metallic  grains  in  the 
cermet  is  obtained  and  the  result  is  shown  in  Figure  1. 

As  seen  from  the  figure  the  effect  of  taking  the  average  is 
to  smooth  out  the  quantum  structure  of  a  single  particle  when 
using  the  QSE  approximation. 
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FIGURE  1  Unaveraged  and  averaged  quantum  dielectric 
constant  of  Ag  metallic  grains  for  a  grain  distribution 
with  average  radius  of  30  A  and  a  geometrical  standard 
deviation  of  1.3 


FIGURES  2a  and  2b  Experimentally  obtained  real  and  imagi¬ 
nary  parts  of  the  index  of  refraction  for  Au/A1203  (90% 
A1203)  cermet  compared  with  BD  and  BQ  models  using  20  A 
average  radius  and  a  limited  mean  free  path  of  10  A. 
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Small  metallic  inclusions  in  cermets  are  known  to  grow  by 
coalescence  of  smaller  inclusions  which  difuse  on  the  substrate, 
as  a  result,  a  significant  number  of  imperfections  are  trapped 
within  the  metallic  grain  during  the  formation  process.  The 
effect  of  these  imperfections  is  considered  introducing  a  restric¬ 
ted  mean  free  path  which  is  smaller  than  the  bulk  mean  free  path 
when  applying  the  QSE  calculation. 

We  now  compare  with  experimental  data  the  results  obtained 
when  the  QSE  dielectric  constant  is  introduced  into  the  calcula¬ 
tion  of  the  effective  dielectric  constant  of  the  cermet.  The 

theoretical  calculations  will  be  compared  with  the  measured 

2 

refractive  index  for  the  Au/A1203  cermet  obtained  by  Craighead  . 
This  experimental  data  was  selected  because  it  is  made  on  cermet 
with  metallic  grains  of  average  radiuses  less  than  SO  A,  and 
consequently,  the  QSE  method  should  contribute  significantly  to 
the  complex  refractive  index. 

Figures  2a  and  2b  compare  the  experimentally  measured  complex 
refraction  index  for  Au/A^O^  cermet,  with  the  Bruggeman's  theory 
using  Drude's  model  (BD)  and  Bruggeman's  theory  using  the  QSE 
correction  (BQ) .  The  average  radius  used  for  the  calculations  was 
20  A,  and  the  value  for  the  limited  mean  free  path  was  10  A.  As 
can  be  observed  from  the  figure  the  BD  model  gives  good  results 
for  the  real  part  of  the  index  of  refraction  but  exaggerates  the 
magnitude  of  the  observed  resonance  in  the  imaginary  part  of  the 
index  of  refraction  (K).  The  BQ  model  greatly  improves  on  the  BD 
for  K  by  reducing  the  magnitude  of  the  peak  and  approximating  much 
better  the  actual  values  of  K  in  the  infrared  region  of  the  spec¬ 
trum.  The  BQ  model  gives  a  smaller  value  for  the  real  part  of  the 
index  of  refraction  N  than  the  BD  model. 

CONCLUSIONS 

-  ■  t- 

The  classical  approach  to  the  theory  of  the  optical  properties  of 
cermets  has  traditionally  used  Drude's  model  to  obtain  the  complex 
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dielectric  constant  of  the  metallic  component  corrected  by  an 
effective  relaxation  time  that  taken  into  account  surface  scatter¬ 
ing,  but  fails  to  consider  quantum  effects.  This  approach  is 
found  to  work  when  the  radius  of  the  metallic  inclusions  are 
larger  than  100  A.  We  have  shown  that  when  the  radiuses  of  the 
metallic  grains  are  smaller  than  a  100  A  QSE  become  an  important 
correction  to  the  effective  dielectric  constant  of  the  metallic 
grains,  when  calculating  the  optical  properties  of  cermets.  QSE 
predicts  discrete  absorbtion  peaks  (see  Part  I)  that  are  not 
observed  experimentally  in  the  optical  measurements  of  cermets. 
Micrograph  of  dielectric  rich  cermets  shows  that  the  metallic 
inclusions  obey  a  log-normal  size  distribution.  When  the  results 
for  the  QSE  for  metallic  grains  are  averaged  over  this  type  of 
distribution  the  absorbtions  are  smoothed  out  yielding  an  average 
single  particle  with  a  smooth  complex  dielectric  constant,  thus, 
explaining  why  no  structure  is  observed  experimentally. 

The  ability  to  reproduce  the  experimentally  observed  optical 
properties  of  cermets  are  considerably  improved  when  the  dielec¬ 
tric  constant  of  the  metallic  constituent  is  corrected  for  Quantum 
Size  Effects. 

We  will  develop  a  multiple  scattering  correction  to  the 

3 

effective  dielectric  constant  of  the  cermet  in  a  following  paper. 
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Abstract  Improvements  to  the  theory  of  the  optical  pro¬ 
perties  of  real  cermets  are  made  by  performing  multiple 
scattering  corrections  to  the  effective  medium  dielectric 
constant.  A  multiple  scattering  model  that  considers 
correlation  effects  between  metallic  grains  in  a  systematic 
way  is  discussed.  The  results  using  both  mean  field  and 
multiple  scattering  theories  are  compared  against  experimen¬ 
tal  data  for  Au/A^O^  and  Ag/MgO  cermets. 


INTRODUCTION 


Search  for  adequate  theoretical  models  that  will  effectively 

predict  the  optical  properties  of  composite  materials  is  a  problem 

1  2 

that  has  been  studied  for  many  years  ’  but  no  satisfactory  solu¬ 
tion  to  the  problem  has  yet  been  found. 

Cermets  are  special  types  of  composite  materials  made  of  two 
constituents,  one  metallic  and  the  other  dielectric.  The  charac¬ 
teristic  cermet  possesses  a  granular  structure  where  the  typical 
diameter  of  the  grains  is  of  the  order  of  a  100  A.  Given  the 
small  size  of  the  microstructure  of  the  typical  cermet  it  is 
always  possible  to  describe  its  infrared  and  visible  optical 
properties  in  terms  of  an  effective  medium  dielectric  constant. 

Mean  field  theories  have  been  developed  that  express  the 
effective  dielectric  constant  of  the  cermet  in  terms  of  the  known 
parameters  of  the  constituent,  to  wit,  the  radius  of  the  grains, 
relative  volume  concentration,  and  dielectric  constant  of  the  two 

constituents.  The  best  known  and  most  basic  effective  medium 

1  2 

theories  are  the  Maxwell-Garnett  -  and  the  Bruggeman  theories. 
These  theories  assume  the  electrostatic  field  approximation 
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spherical  shapes  or  other  similar  geometric  structure  for  the 

grains  and  a  completely  random  microstructure.  Many  modifications 

of  these  theories  have  been  developed,  correcting  for  different 
3  4 

shape  factors  ,  island  size  distributions  ,  and  effective  depola¬ 
rization  factors^.  Some  attemps  have  been  made  to  include  in  an 
indirect  way  correlation  effects  between  the  grains  in  the  micro¬ 
structure  by  considering  coating  effects**  or,  in  some  few  cases, 
by  introducing  a  correlation  function  for  the  microstructure7. 

The  objective  of  this  work  is  to  briefly  discuss  a  method 
that  explicitly  takes  into  account  correlation  between  metallic 
grains  and  corrects  for  multiple  scattering  effects  in  calculating 
the  effective  dielectric  constant  of  the  cermet. 

Multiple  scattering  correction  will  be  compared  with 
Bruggemann's  model  against  experimental  data  for  a  Ag/MgO.  For 
AU/AI2O3  cermet  the  metallic  grain  sizes  are  significantly  smaller 
than  100  A  and  the  model  shall  be  used  together  with  quantum  size 

g 

effects  and  compared  with  mean  field  theories  and  experimental 
data  to  demonstrate  that  these  corrections  provide  a  better  fit 
with  the  results  obtained  from  optical  measurements  of  real 
cermets. 


Multiple  Scattering  Correction 


Standard  method  for  growing  cermet  films  such  as  sputtering  and 
co-evaporation  produces  a  film  growth  process  dominated  by  dif¬ 
fusion  and  coalescence  phenomenas.  Due  to  these  mechanisms  the 
positions  of  adjacent  grains  are  strongly  correlated  with  a 
gradual  decrease  in  the  correlation  between  grains  as  the  distance 
between  them  increases.  As  a  first  approximation  for  grain-grain 
correlation  an  exponential  type  of  correlation  function  can  be 
utilized.  A  more  careful  development  of  a  theory  for  correlation 

effects  in  cermets  worked  out  by  our  group  corroborates  that  the 

9 

exponential  approximation  is  an  excellent  one  .  The  validity  of 
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octal  grains  are  totally  uncorrclated .  Since  this  is  not  the  case 
for  real  cermets,  we  have  removed  this  restriction  and  constructed 
a  multiple  scattering  theoretical  model  that  effectively  takes 
into  account  these  correlations . 

The  fluctuations  in  the  dielectric  constant  due  to  the 
microsctructure  are  taken  to  be  deviations  from  an  average 
dielectric  constant  that  we  will  take  to  be  Bruggeman's  mean  field 
result.  That  is 

e(r,u>)  =  e0(u>)  ♦  6e(r,u»)  (1) 

,  where  Eo(ui)  is  the  Bruggeman's  value  for  the  dielectric  constant 
of  the  cermet  and  6e(r,ui)  is  the  fluctuation  due  to  the 
microstructure. 

Our  theory  uses  a  perturbation  expansion  of  the  propagating 
electric  field  in  the  cermet  expressed  as  a  power  series  of  the 
fluctuating  part  of  the  dielectric  constant. 

Starting  from  the  electric  field  equation 

V  x  V  x  £  -  (u>/c)2peE  -  m’1Vp  *  V  *  l  -  0  (2) 

This  equation  should  be  approximated  and  expressed  in  terms 
of  the  fluctuating  part  of  the  dielectric  constant  as: 

V  x  V  x  j?  -  k022  =  Kg  (6 e/e0)  E  (3) 

,  where  Kq  is  the  effective  propagation  constant  in  the 
Bruggeman's  medium.  Then  the  field  equation  may  be  written  as  an 
integral  equation.  Using  the  operator  form, 

G  =  Gc  +  G0VG  (4) 

,  where  G  and  Gq  are  the  Green's  tensors  associated  with  the  field 
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£  and  Che  field  in  Che  Bruggeman  medium  respeccively .  V  is  Che 
perCurbation  Cera  due  Co  the  flucCuaCions  fie. 

An  average  is  Chen  perforaed  Co  obcain  a  macroscopic  field 
equaCion  keeping  only  terms  involving  pair  correlation  functions 
oi  the  fluctuating  part  of  Che  dielectric  constant  and  using  the 
so  called  First  Smoothing  Approximation  (FSA)1  ' 11 . 

<G>=G0  ♦  G0  <VG0V>  <G>  (5) 

The  equation  obtained  is  an  infinite  series  corresponding  only  to 
scattering  processes  where  two  point  correlations  are  involved. 
Diagramatically 

<G>  =  -  +  +■  •  •  •  (6) 

,  where  symbol  -  represents  G0,  •  is  V,  and  the  diagram 
contains  the  pair  correlation  function  <fie(ri)fie(r2)>. 

12 

Following  the  procedures  used  by  F.  Karal  and  J.  Keller  a 
plane  wave  solution  is  then  assumed  for  the  macroscopic  field  in 
the  equation  resulting  from  the  FSA  and  the  effective  dielectric 
constant  corrected  for  multiple  scattering  effect  is  then 
obtained.  To  get  the  final  result,  the  pair  correlation  function 
was  chosen  to  be  an  exponential  decay  type  as  we  pointed  out 
above ,  and 

<fie(?l)fie(?2)>  =  A  e"r/L  (7) 

,  where  L  is  the  correlation  length,  r=  r2  -  rt  ,  and  A  was  chosen 
as  the  square  of  the  standard  deviation  of  the  fluctuations. 

The  final  result  expresses  the  effective  dielectric  constant 
of  the  composite  in  terms  of  the  dielectric  constant  obtained 
using  Bruggeman  model  and  the  correlation  length  L.  The  region  of 
interest  for  the  optical  properties  of  cermets  lies  between  0.5ev 
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and  4.0ev  in  the  electromagnetic  spectrum,  in  this  range  the  final 
result  should  be  expressed  as: 


£ 


2_  1+ia  2a2 

3  a2  l-i2a 


,  where  a2  =  £oiu2I2/c2. 


(8) 


Application  of  Multiple  Scattering  Corrections  to  Real  Cermets 

We  will  now  apply  the  developed  multiple  scattering  corrections 

theory  (M)  to  perform  the  calculation  of  the  optical  properties  of 

real  cermets  and  compare  the  results  with  those  obtained  using 

Bruggeman' s  theory  (B) .  The  predictions  of  both  theoretical 

models  B  and  M  will  be  compared  against  experiisental  data  obtained 
13 

by  H.  Craighead  for  the  real  (N)  and  imaginary  part  (K)  of  the 
index  of  refraction. 

Figures  la  and  lb  show  the  experimental  results  for  Ag/MgO 
cermet  with  a  concentration  of  80%  of  the  dielectric  constituent 
compared  with  B  and  M  models.  The  dielectric  constant  of  the 
metallic  grains  is  an  important  paraawter  included  in  the  B  and  M 
models.  In  that  cermet  we  have  calculated  it  using  Drude's  model 
corrected  by  a  restricted  mean  free  path  due  to  impurities  and 
defects  into  the  metallic  grains.  Interband  contributions  are 
also  taken  into  account.  Another  important  parameter  in  the  M 
model  is  the  correlation  length  that  was  125  A  in  that  calculation 

As  can  be  observed  from  the  figures,  M  obtained  using  N  model 
reproduces  B  model  when  \<lpa  but  predicts  lower  values  for  larger 
wavelengths  that  are  also  obtained  with  the  experimental  data.  B 
model  predicts  higher  values  for  the  imaginary  part  of  the  index 
of  refraction  K  than  those  obtained  experimentally  in  all  the 
region  of  interest  but  M  model  performs  a  very  good  fit  with  the 
experimental  data  in  almost  all  the  region.  Only  a  slight  differ¬ 
ence  is  found  around  0.6pm.  In  Figure  1  the  best  fit,  to  date,  is 
found  between  experiment  and  theory  for  this  cermet. 
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FIGURE  la  and  lb  Experimental  real  and  imaginary  parts  of 
the  index  of  refraction  for  Ag/MgO  (80%  MgO)  cermet  compared 
with  B  and  M  models  using  a  restricted  mean  free  path  of  5  A 
and  a  correlation  length  of  125  A. 


FIGURE  2a  and  2b  Experimental  real  and  imaginary  parts  of 
the  index  of  refraction  for  AU/AI2O3  (84%  A1203)  cermet 
compared  with  BO,  BQ  and  MQ  with  80  A  correlation  length,  30 
A  average  radius  and  limited  mean  free  path  of  10  A. 
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Figures  2a  and  2b  compare  the  theoretical  models  with 

experimental  results  for  the  Au/A1203  cermet  with  84%  volume 

fraction  of  the  dielectric  constituent.  In  this  cermet  the 

metallic  grains  are  smaller  than  100  A  and  quantum  corrections  to 

8 

the  dielectric  constant  of  the  metallic  grains  become  important  . 
To  obtain  the  best  results  with  the  multiple  scattering  model 
using  quantum  size  corrections  (MQ) ,  the  restricted  mean  free  path 
was  taken  to  be  10  A,  and  the  correlation  length  to  be  80  A.  In 
the  figures  we  compare  Bruggeman  model  using  Drude's  dielectric 
constant  with  interband  contributions  for  the  metallic  grains 
(BD),  Bruggeman  with  quantum  size  corrections  (BQ)  and  MQ  against 
experimental  data.  The  BD  model  again  predicts  a  higher  peak  in  K 
than  the  one  observed  experimentally  and  the  values  obtained  for 
this  model  are  always  larger  than  the  experimentally  observed  ones 
for  both  N  and  K.  The  BQ  results  improved  on  the  BD  model  fitting 
much  better  the  experimental  N  and  K.  The  MQ  model  greatly 
improves  on  the  BD  and  BQ  models  by  narrowing  and  decreasing  the 
peak  for  K  and  giving  better  agreement  with  the  experimental 
result  in  the  region  0.8y m  <  A  <  2.5|Jm  while  coinciding  with  the 
BQ  model  for  smaller  wavelength. 

CONCLUSION 

The  processes  that  involve  the  cermet  formation  make  that  the 
position  of  adjacent  grains  are  strongly  correlated.  A  multiple 
scattering  theory  that  takes  into  account  these  correlation 
effects  in  a  systematic  way  is  needed  to  predict  the  optical 
properties  of  cermets. 

A  multiple  scattering  expression  for  the  effective  dielectric 
constant  of  a  cermet  is  developed  as  a  correction  to  the  Bruggeman 
mean  field  theory  and  used  to  predict  the  values  of  the  effective 
index  of  refraction  of  real  cermets. 

Comparison  against  the  experimental  data  shows  that  the 
theory  considerably  improves  the  Bruggeman' s  model.  When  the 
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sizes  of  the  metallic  grains  of  the  cermet  are  smaller  than  100  A, 
quantum  size  corrections  become  important  and  it  must  be 
considered  together  with  scattering  corrections,  then  the  best  fit 
between  experiment  and  theory  is  obtained. 

Multiple  scattering  effects  are  important  when  the  volume 
fraction  of  the  metallic  constituent  in  the  cermet  is  greater  than 
10%  and  grain-grain  correlations  become  significant.  Work  is  in 
progress  to  further  improve  the  multiple  scattering  theory  by 
taking  into  account  renormalization  effects  within  the  T-Matrix 
formalism. 
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multiple  scattering  renormalized  t-matrix  theory  to  the  oielectric  constant  of 

NON'HOMOGENEOUS  THIN  FILMS.  Manual  Gomez  and  Luis  Fonseca.  Physics 
Oepartment,  University  of  Puerto  Rico,  Rio  Piedras,  Puerto  Rico. 

I.  INTRODUCTION 

The  inhomogeneous  media  of  granular  metal  films  have  been  extensively 

1 

studied  because  of  their  importance  to  technological  applications  ‘  The 

optical  property  of  films  composed  of  random  arrays  of  metallic  grains 

embedded  in  a  dielectric  matrix  (cermet)  have  been  studied  for  their  spectral 

2 

selectivity  which  make  them  promising  selective  surfaces  . 

In  the  range  of  the  visible  and  near  infrared  regions  of  the  spectrum 
where  the  wave  length  of  the  incident  light  is  large  compared  to  the  radiuses 
of  the  metallic  grains,  the  effect  of  the  microstructure  is  averaged  and  the 
film  behaves  as  an  effective  homogeneous  medium.  Several  theoretical  models 
have  been  proposed  to  obtain  the  effective  dielectric  constant  of  cermets  in 
terms  of  the  dielectric  constant  of  the  constituents  and  parameters  describing 
the  microstructure,  such  as  island  radiuses,  relative  concentrations  of  the 
constituents,  and  the  shape  of  the  islands.  Two  basic  and  widely  used  models 

3 

are  the  Maxwel  1 -Garnett  (MG)  and  the  Bruggeman  (B)  models  .  The  MG  model  is 
asymmetric,  since  it  considers  the  cemtet  as  composed  of  spheres  of  the  less 
abundant  constituent  embedded  in  a  matrix  of  the  most  abundant  material.  The 
effective  dielectric  constant  is  then  obtained  by  performing  a  volume  average 
over  the  local  fields.  While  the  B  model  is  symmetric  since  it  considers  the 
two  constituents  materials  as  spheres  embedded  in  an  effective  medium,  a 
volume  average  is  also  performed  to  obtain  a  self-consistent  expression  for 
the  effective  dielectric  constant.  In  both  theories  a  resonance  in  the 
spectral  response  of  the  cermet  is  obtained,  but  neither  agrees  with  the 
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experimentally  measured  complex  dielectric  constant.  For  example,  comparison 
between  the-  experimentally  obtained  imaginary  part  of  the  effective  index  of 
the  refraction  and  the  MG  model  reveals  that  the  theoretical  resonance  is  too 
high  and  narrow  and  in  many  cases  blue  shifted  with  respect  to  experimental 
results.  Several  papers,  like  those  of  Granqvist*  show  that  the  magnitude  and 
position  of  the  resonance  is  a  strong  function  of  the  inclusions'  shape. 
Corrections  to  the  dielectric  constant  of  metallic  grains  for  size  effects 
have  been  extensively  used  to  reduce  the  magnitude  of  the  resonance,  but  they 
fail  to  produce  enough  red  shift  and  broadening  to  explain  experimental  data. 

II.  THE  THEORY 

We  consider  the  cermet  as  a  random  array  of  metallic  grains  embedded  in  a 
dielectric  matrix.  The  spatial  fluctuation  of  the  dielectric  constant  can  be 
envisioned  as  consisting  of  a  constant  dielectric  constant  eQ  (e.g.  the 
dielectric  constant  of  the  more  abundant  material  or  a  judiciously  select  mean 
value  dielectric  constant)  and  a  stochastic  term  5c  that  accounts  for  the 
spatial  fluctuations  of  the  dielectric  constant,  i.e. 

c(r)  =  eQ  *  6e(r) 

With  this  expression  the  field  equations  for  the  propagating  electromagnetic 

field  can  be  separated  into  a  non-stochastic  and  a  stochastic  fluctuating  term. 

The  fluctuating  tens  will  be  considered  as  the  perturbing  potential 

responsible  for  the  scattering  processes.  The  resulting  field  equations 

cannot  be  solved  exactly  and  approximation  methods  must  be  used. 

s  e 

The  First  Smoothing  Approximation  (FSA)  ’  considers  only  the 
contribution  of  two  point  correlations  to  the  multiply  scattered  field.  This 
approximation  has  previously  been  applied  to  obtain  the  effective  dielectric 

7  * 

constant  of  a  general  random  medium  and  more  recently  to  cermets  .  Our 
recent  work  for  cermets  showed  that  the  FSA  approximation  corrects  the 
magnitude  of  the  resonance  but  fails  to  predict  a  significant  spectral  shift, 


suggesting  that  renormalization  corrections  to  the  scattering  terms  must  be 

included  in  order  to  obtain  better  agreement  with  experiments.  In  order  to 

e 

apply  renormal ization  techniques  ,  the  T-matrix  multiple  scattering  of  an 
individual  grain,  is  redefined  as  the  new  renormalized  potential  symbolically 
represented  by  $  ,  and  responsible  for  scattering  processes  in  the  medium.  A 
new  Green's  function  is  then  generated  by  taking  into  account  multiple 
scattering  effects  from  these  new  single  scatterer  potentials  in  the  absence 
of  correlation.  These  Green's  function  associated  with  <£>  and  represented 

diagramatical ly  by  -  then  defines  an  effective  medium,  that  because  it 

is  uncorrelated,  we  will  call  the  zeroth  order  effective  medium  characterized 
by  a  dielectric  constant  eQ.  In  terms  of  the  new  renormalized  Green's 
function  and  potential  an  averaged  field  equation,  that  takes  into  account 
pair  correlation  functions,  is  written  to  yield  the  Renormalized  First 
Smoothing  Approximation  (RFSA), 


or,0*or  or  3 


where  <?Q>  is  the  statistically  averaged  field  propagating  in  the  uncorrelated 
medium  and  <£>  the  effective  total  field.  The  second  term  in  the  r.h.s.  part 
of  this  self-consistent  equation  represents  scattering  processes  where  the 
effective  total  field  <£>  is  scattered  by  pairs  of  correlated  islands  or  and  3 
located  respectively  at  r.  and  tn  order  to  apply  this  formalism  an 

explicit  mathematical  representation  of  the  diagram  in  the  second  term  of  the 

9 

r.h.s.  of  eqn.  (1)  must  be  obtained.  Peterson  and  Strom  developed  a 
formalism  to  obtain  the  T-matrix  from  an  array  of  n-scattering  objects  with 

well  defined  spatial  coordinates.  Their  formalism  is  an  extension  to  multiple 

10 

number  of  scatterers  of  that  developed  by  Waterman  for  a  single  scatterer. 
In  Waterman's  formalism  the  T-matrix  for  a  single  scatterer  is  obtained 


4 


expanding  all  tha  local  fields  in  terms  of  Spherical  Harmonics,  and  Bessel  or 

Hankel  functions.  The  term  describing  this  double  scattering  process  in 

Peterson  and  Strom's  formalism  is  given  by: 

®— ®  -  R(fj)  T®  TP  RCPj.-f.)  R(-?j)  ,  (2) 

where  T® ,  TP  are  the  single  scatterer  T-matrices  of  islands  a  and  0, 

respectively,  expressed  with  their  center  in  the  origin  of  the  coordinate 

system.  The  formalism  then  translates  the  particle  through  matrices  R  and  a 

to  their  correct  position  r .  and  r..  Expressions  for  T®  are  described  in 
9  to 

several  articles  ’  and  expressions  for  matrices  R  and  a  can  be  found  in 
Peterson  and  Strom's  article.  Using  eqn.  (2)  the  field  equation  can  now  be 
expressed  as: 


<£(?)>  = 


<Eorr)> 


1/V 


I 

a,0*cr 


// 


g(r 


i'V 


T®  oCr.-fj) 


(3) 


x  TP  R(r‘j-ri)  R(-Fj)  <t>  dr*,  d^ 

where  g  (r^.fj)  is  the  two-point  correlation  function  and  V  is  the  volume  of 
the  film.  Considering  the  spheres  as  impenetrable  and  the  medium  as  isotropic 
exponential  correlation  functions  can  be  defined  as: 

3  0(r-2a)  exp  [- |r.  -  r.|/  L] 

where  0  is  the  Heaviside  function  and  L  the  correlation  length.  To  simplify 

the  calculation  we  assume  all  the  metallic  grains  to  be  spherical  in  shape 

and  of  equal  radiuses  which  allow  us  to  write  T®=TP  and  make  them  diagonal 

2  2 

matrices.  Assuming  a  microscopic  volume  we  can  change  (1/V  )I  by  5  ,  where  6 

Is  the  density  of  Islands  In  the  cermet.  The  field  <?>  and  <EQ>  will  be 

expanded  In  terms  of  the  basis  functions  {Re(jin(Kor)}  that  are  solutions  to  the 

io  - 

Helmholtz  equation  described  by  Waterman  The  ^  are  known  as  the 

n 

elementary  fields  and  KQ  is  the  propagation  constant  of  the  electromagnetic 


5 


wave  In  the  zeroth  order  effective  medium.  In  terms  of  this  base  the  fields 
can  be  expressed  as: 

<fo>sIdnK(l(of)  and  <?<?)>  =ZXn  Re*n(Ko?)  • 

n  n 

While  the  total  effective  field  <E>  in  the  second  term  of  the  r.h.s.  is 
expanded  In  a  base  centered  at  one  of  the  correlated  islands: 

<i>  =£en 

5* 

Substituting  these  expansions  for  the  fields  in  eqn.  (3)  and  performing  the 
integration  over  the  relative  coordinate  between  the  correlated  islands  we 
obtained  : 


Z*n  R«n<Kof>  ’Zdn  *  5TMn  /  8n  af3  ■  <4> 

a  a  a  ' 


Matrix  elements  M  contains  all  the  correlations  effects  and  can  be  written  in 
n 

terms  of  the  elements  of  the  single  scatterer  T-matrices.  Since  the  wave 
length  of  Incident  radiation  relevant  to  our  problem  is  A.  >  .3  pm  and  the 
radiuses  of  Islands  for  typical  cermets  are  of  the  order  of  .01  pm,  the  long 
wave  limit  can  be  utilized,  thus  reducing  the  expansion  of  the  electric  fields 
to  the  first  few  terms.  In  this  long  wave  limit  only  the  first  six  terms  of 
these  matrices  are  necessary,  that  correspond  to  a  dipolar  approximation. 

To  perform  the  integral  over  the  coordinate  we  must  use  the  properties 
of  the  translation  matrices  cj11  to  translate  the  $n  from  position  r^  to  the 
origin, 

VV^V3’IVn<W  •  C5> 

h* 


and  all  terms  In  eqn.  (4)  are  now  written  in  the  same  orthogonal  base  and  is 


equivalent  to  the  following  set  of  scalar  equation: 

X  ■  d  *  f  8  ,  ann,  dr  . 
n  n  I  n  n  nn  j 


tv  J 
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Since  we  are  in  the  long  wave  limit  where  the  electromagnetic  wave  tends  to 
average  the  microstructure  at  the  scale  of  island  sizes  we  can  assume  that  the 
wave  propagates  in  the  effective  medium  as  a  plane  wave.  The  coefficient  d 

n 

can  be  written  in  terms  of  KQ,  while  the  coefficient  representing  <?>  can  also 
be  expressed  in  terms  of  the  effective  propagation  constant  K  of  the 
correlated  medium  as: 

dn  *  exPt,Ko‘?1  and  Xn  =  XJ  exPf’*,?]  and  Bn  =  X*  exp(ii?-(r-r .)] 

J 

In  terms  of  these  expressions  the  field  equation  now  becomes: 

X*  .xp[t*.p]  >  d”  «xp[1Ko.?]  .  4*  £«„,  Xn,  f  ,,p[(R.(r-f.)]  onn,(K/  )  dr 

n'  J  .  ,  J 


•m  -• 

r‘rjl 


Following  the  procedure  described  by  Varadan  et  al.  we  can  use  the  scalar 
Helmholtz  operator  in  the  zeroth  order  approximation  to  simplify  the  above 
equation  and  reduce  the  volume  integral  to  a  surface  integral: 

<  *  I  <•  'op'  '  (6) 


where 

2  2  2 
Xnn‘  =  6  /(K  'Ko  J 


f  {a nnl(‘Ko? 
I  r|  >  a 


)  3exp(iK*r)/3r  -  exp(iK*r)  3ann, (-KQr)/3r}  ds 


Eqn.  (6)  represents  a  system  of  homogeneous  coupled  equations  with  unknown 
coefficients.  Solving  the  associated  secular  equation  we  obtained  a  new 
dispersion  relationship  that  takes  into  account  pair  correlations  in  the 
scattering  medium.  In  terms  of  the  dielectric  constants  this  expression  can 


be  expressed  as: 


5  *  l0  (1  ♦  2Fk>E0)/(l-FkSf0) 


where  e  is  the  new  effective  dielectric  constant  in  the  correlated  medium  of 
the  film  obtained  from  the  RFSA  model,  k  =  w/c  is  the  propagation  constant  of 
the  electromagnetic  wave  in  vacuum,  and  F  contains  all  the  effects  due  to 
correlations  between  islands  and  can  be  expressed  as: 

F=  22/50  f*  exp(-2a/L)  (l/*2al)  ((ep-l)/(er«-2)]* 
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wh«r*  (a)  is  the  radius  of  the  inclusions,  (f)  the  relative  concentration  of 
metal,  and  er  is  the  ratio  of  the  dielectric  constants  of  the  metallic  islands 
and  the  external  medium. 

III.  RESULTS  AND  DISCUSSION 

The  result  obtained  in  eqn.  (7)  permits  the  calculation  of  the  complex 

effective  dielectric  constant  of  cermets  taking  into  account  correlation 

effects,  in  terms  of  basic  physical  parameters  obtainable  from  experimental 

data.  When  a  metal  concentration  is  small  the  correlation  length  tends  to 

vanish  and  we  recover  from  eqn.  (7)  the  zeroth  order  approximation  eQ,  the 

effective  dielectric  constant  for  the  uncorrelated  effective  medium.  To  apply 

the  dispersion  relationship  the  MG  model  is  used  for  the  zeroth  order 

dielectric  constant.  Conceptually  MG  model  is  compatible  with  the  T-matrix 

formalism,  in  addition,  Varadan  et  al.,  have  shown  that  in  the  limit  of  no 

12 

correlation  the  T-matrix  calculations  yields  the  MG  result  For  the 

dielectric  constant  of  the  metallic  islands  both  Drude's  model  (D)  with 

relaxation  time  t  corrected  for  size  effect  and  the  Quantum  Size  Effect  model 
1 3 

(QSE)  will  be  used,  both  calculations  are  corrected  for  the  interband 
contributions. 

Fig.  1  shows  the  real  and  the  imaginary  part  of  the  refraction  index  for 
Ag/MgO  80%  MgO  cermet.  In  this  case  the  t  used  in  the  Qrude  model  is  adjusted 
to  give  the  best  fit  to  the  experimental  resonance  by  the  MG  model.  Then  using 
the  same  dielectric  constant  for  the  metal  constituent,  the  RFSA  calculations 
shifts  and  broadens  the  peak  improving  the  MG  calculation  and  giving  a  much 

14 

better  fit  to  experimental  data  obtained  from  H.  Craighead  work  .  Fig.  2 
shows  the  same  experimental  data,  but  now  the  D  model  has  been  substituted  for 
the  QSE  model  in  the  calculation  for  the  metallic  dielectric  constant.  In 
this  case  a  more  realistic  t  gives  the  best  fit  to  experiment  by  the  RFSA 
model  but  the  peak  is  now  narrower.  The  corresponding  mean  free  path  for  the 
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t  used  in  this  figure  is  1QA.  Values  of  up  to  5A  had  to  be  used  by  Craighead 
to  obtain  the  best  fit  with  the  MG  theory  to  his  experimental  data.  We  have 
also  applied  the  model  to  other  cermets  such  as  Ni/Alz  03  obtaining  a  very 
good  fit  with  experimental  data  with  metallic  concentrations  of  up  to  54%. 

IV.  CONCLUSIONS 

We  have  extended  our  FSA  model  calculations,  that  take  into  account 

correlation  effects  among  metallic  islands  in  cermet  films,  to  include 

renormalized  corrections  and  have  shown  that  correlation  effects  are  important 

in  the  broadening  and  shifting  the  resonance  peak  to  obtain  better  fit  with 

a 

experimental  results.  The  results  of  this  and  other  recent  work  leads  us  to 

conclude  that  correlation  effects  and  multiple  scattering  corrections  make 

important  contributions  to  the  optical  properties  of  cermet  materials  and  need 

to  be  included  in  theories  with  predictive  value.  Eqn.  (7)  is  restricted  to 

spherical  shapes  and  long  wave  limit,  but  the  formalism,  as  presented  here, 

can  be  extended  to  include  other  island  shapes  and  to  permit  corrections  for 

higher  multipolar  contributions  to  the  field  equation  that  may  become 

important  when  the  metal  concentration  increases  and  island  proximity  may  make 

these  consideration  important.  The  proposed  dispersion  relation  depends 

strongly  on  the  models  used  for  the  zeroth  order  dielectric  constant  and  used 

for  dielectric  constants  of  the  metallic  grains.  Other  zeroth  order  models 

for  the  effective  medium  dielectric  constant  are  being  studied  as  possible 

starting  points  for  the  application  of  the  dispersion  relationship  presented 

in  this  paper.  These  new  calculations  might  permit  the  use  of  a  longer  t  that 

correspond  to  more  realistic  mean  free  path.  Recently  other  models  have  been 

proposed  for  tQ  that  should  help  moderate  the  height  and  broaden  the  resonance 
is 

with  a  larger  t 
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SUMMARY 

A  multiple  scattering  theory  to  obtain  the  effective  dielectric  constant 
for  a  granular  metal  film  is  developed  considering  correlation  effects.  A 
renormalization  procedure  is  applied  to  obtain  a  field  equation  within  a 
Renormalized  First  Smoothing  Approximation  (RFSA),  using  Tmatrix  formalism. 
From  the  field  equation  a  dispersion  relation  is  obtained  that  is  expressed  in 
terms  of  the  dielectric  constant  of  the  metal  corrected  for  size  effects,  the 
relative  concentration  of  the  constituents,  the  radius  of  the  inclusions  and  a 
correlation  length.  Using  the  Maxwel 1 -Garnett  model  for  the  dielectric 
constant  of  the  uncorrelated  medium,  the  multiple  scattering  theory  shows  that 
correlation  effects  moves  the  resonance  in  the  imaginary  part  of  the  index  of 
refraction  predicted  by  the  Maxwel 1 -Garnett  model  toward  lower  frequencies  and 
broadens  it.  Comparison  of  the  RFSA  with  experiments  shows  better  fit  than 
previous  theories. 
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FIGURE  CAPTION 

Figure  1: 

The  real  part  and  the  imaginary  part  n ^  of  the  refractive  index  as  a 
function  of  the  wavelength  is  drawn  for  an  Ag/MgO  80%  MgO  cermet.  Solid  line 
is  the  result  of  RFSA  calculation,  dashed  lines  that  of  the  MG  model  and  dots 

14 

represent  the  experimental  results  from  H.  Craighead  This  figure  is 

vs 

obtained  using  the  0  model  with  t=  2.14  x  10  sec  for  the  dielectric  constant 
of  the  metallic  grains  in  both  MG  and  RFSA  curves.  The  radius  of  the  islands 
is  taken  as  80A  and  L=6  islands. 

Figure  2: 

The  real  and  imaginary  parts  of  the  refractive  index  are  drawn  for  the 

13 

same  cermet  of  figure  1.  This  figure  is  obtained  using  QSE  with  a 
restricted  mean  free  path  of  10A  for  the  dielectric  constant  of  the  metallic 
grains  in  both  MG  and  RFSA  curves.  The  island  radius  is  125A  and  L=6  islands. 
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Introduction 

Cermets  are  inhomogeneous  materials  consisting  of  immisibles  mixtures  of 
insulators  and  metals.  Extensive  work  has  been  performed  to  develop 
theoretical  models  which  will  correctly  predict  the  optical  properties  of 
cermets  when  the  bulk  dielectric  constant  of  the  constituents  is  known 
experimentally.  When  the  concentration  of  the  dominant  constituent  is  less 
than  98%  the  material  develops  a  distinct  resonance  which  is  not 
characteristic  of  either  one  of  the  constituents.  The  prediction  of  this 
resonance  has  been  the  goal  of  all  theories  developed  for  the  optical 
properties  of  cermets  materials. 

Two  basic  and  widely  used  models  are  the  Maxwel  1  -Garnett  (MG)  and  the 
Bruggeman  (B)  models  .  The  MG  model  is  asymmetric,  since  it  considers  the 
cermet  as  composed  of  spheres  of  the  less  abundant  constituent  embedded  in  a 
matrix  of  the  most  abundant  material.  The  effective  dielectric  constant  is 
then  obtained  by  performing  a  volume  average  over  the  local  fields.  While  the 
B  model  is  symmetric,  since  it  considers  the  two  constituent  materials  as 
spheres  embedded  in  an  effective  medium,  a  volume  average  is  also  performed  to 
obtain  a  self-consistent  expression  for  the  effective  dielectric  constant.  In 
both  theories  a  resonance  in  the  spectral  response  of  the  cermet  is  obtained, 
but  neither  agrees  with  the  experimentally  measured  complex  dielectric 
constant.  For  example,  comparison  between  the  experimentally  obtained 
imaginary  part  of  the  effective  index  of  refraction  and  the  MG  model  reveals 
that  the  theoretical  resonance  is  too  narrow  and  in  many  cases  blue  shifted 
with  respect  to  experimental  results. 

The  objective  of  this  paper  will  be  to  show  that  corrections  due  to 
correlation  between  the  position  of  the  metallic  islands,  resulting  from  the 
formation  process  of  the  cermet,  can  move  and  broaden  resonance  obtained 
by  the  traditional  mean  field  theories  of  Maxwell -Garnett  and  Bruggeman.  In 
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order  to  perform  this  correction  a  multiple  scattering  theory  has  been 
developed  that  starts  from  an  effective  dielectric  constant  that  is  calculated 
assuming  no  correlation  between  the  metallic  islands  and  labeled  in  this  work 
as  the  zeroth  order  approximation.  Non-renormal ized  and  renormalized  theories 
will  be  discussed  in  which  correlation  effects  are  taken  into  account.  Also 
corrections  for  size  effects  on  the  dielectric  constant  of  the  metallic  grains 
will  be  studied. 

Quantum  Size  Effects 

In  order  to  calculate  the  effective  dielectric  constant  of  cermets  it  is 

necessary,  in  general,  to  correct  the  bulk  dielectric  constant  of  the  minority 

constituent  for  size  effects.  Normally  the  insulator  has  a  dielectric 

constant  that  is  essentially  non- frequency  dependent  in  the  spectral  range  of 

interest,  while  the  metallic  constituent  has  a  strong  frequency  dependence  in 

the  same  region.  When  the  metal  is  the  minority  constituent,  its  bulk 

dielectric  constant  must  be  corrected  for  size  effects.  Traditionally  Drude's 

model  together  with  experimentally  obtained  interband  contributions  is  used  to 

calculate  the  dielectric  constant  of  the  metal.  In  this  procedure,  the  mean 

free  path  of  the  electrons  in  the  metallic  island  is  corrected  for  boundary 

scattering  at  the  interface  between  the  island  and  the  external  medium.  Many 
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researchers  have  observed  that  in  order  to  obtain  a  reasonably  good  fit 
between  theory  and  experiment,  it  is  necessary  to  use  mean  free  paths  that  are 
considerably  smaller  than  the  ones  attributed  only  to  the  size  of  the  island. 
These  corrections  could  be  explained  as  the  combination  of  boundary  scattering 
at  the  metallic  grains  surface  and  scattering  inside  the  island  due  to 
internal  imperfections.  The  Drude  theory  with  these  corrections  is 
satisfactory  for  metallic  grains  of  more  than  100A.  However,  when  the  size  of 
the  metallic  particles  become  comparable  to  the  electron  wavelength,  quantum 
corrections  become  mandatory  if  a  realistic  dielectric  constant  is  going  to  be 
obtained. 

To  solve  this  problem  several  works  have  used  the  Quantum  Box  Method 
within  the  Random  Phase  Approximation  to  calculate  the  dielectric  constant  of 
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the  metallic  particles.  Kawabata  and  Kubo  have  corrected  the  dielectric 
constant  for  size  effects  assuming  that  the  spectrum  of  the  quantum  levels 

4 

form  essentially  a  continuum,  while  Wood  and  Ashcroft  carefully  reworked  the 
Quantum  Size  Method,  fully  taking  into  account  the  spectrum  resulting  from  the 
box  quantization  and  obtained  a  dielectric  constant  for  the  particles. 
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Recently  this  formalism  has  been  applied  by  the  authors  of  this  article 
to  real  cermets  in^the  region  where  the  metallic  inclusions  in  the  cermets  are 
smaller  than  100A  and  better  results  have  been  obtained  for  the  dielectric 
constant  than  when  Orude's  model  is  used. 

Figure  1  compares  the  imaginary  part  of  the  dielectric  constant  taking 
into  account  quantum  size  effects  following  Wood  and  Ashcroft  (QSE)  and 
Kawabata  and  Kubo  model  (KK)  with  the  one  calculated  using  Orude's  model.  In 
the  calculation  of  the  Orude's  model  the  size  of  the  silver  particles  are 
taken  to  be  30A  and  no  interband  contributions  to  the  dielectric  constant  have 
been  included. 

The  multiple  peaks  appearing  in  the  QSE  calculations  are  due  to  the 
discrete  energy  spectrum  of  an  electron  in  a  box.  Experimentally  the 
resonances  observed  in  the  dielectric  constant  of  an  isolated  metal  grain  due 
to  quantum  size  effect  are  not  observed  in  the  effective  dielectric  constant 
of  a  real  cermet.  The  disappearance  of  these  peaks  is  due  to  the  fact  that 
the  metallic  grains  in  a  real  cermet  do  not  have  a  unique  size,  but  instead 
have  a  distribution  of  sizes  as  can  be  verified  by  micrographs  of  the 
materials.  For  this  reason  the  dielectric  constant  obtained  from  the  QSE 
calculations  for  one  isolated  particle  were  averaged  over  a  distribution  that 
realistically  represented  variation  in  size  in  actual  cermets.  Averaging  the 
dielectric  constant  of  one  metallic  grain  with  a  log-normal  distribution 
function,  an  effective  dielectric  constant  for  the  metallic  grains  in  the 
cermet  is  obtained.  This  average  QSE  result  is  then  used  to  calculate  the 
effective  medium  dielectric  constant  using  the  B  and  MG  mean  field  theories, 
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Figure  1.  Comparison  between  the  Imaginary  part  of  the  dielectric 
constant  as  function  of  frequency  of  the  QSE.O,  end  kk 
models. 
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and  multiple  scattering  corrections  theories.  To  take  into  account  the  effect 
of  the  imperfections  we  will  introduce  a  restricted  mean  free  path  smaller 
than  the  bulk  mean  free  path  when  applying  the  QSE  calculation  to  obtain  an 
effective  complex  dielectric  constant. 

Multiple  Scattering  Theory:  The  First  Smoothing  Approximation  (FSA) 

When  thin  films  of  cermets  are  grown  by  cosputtering  or  coevaporation 
methods,  nucleation  processes  are  responsible  for  the  formation  of  metallic  or 
dielectric  islands  in  the  matrix  of  the  correspondingly  more  abundant 
constituent.  Since  the  nucleation  process  is  governed  by  the  diffusion 
effects,  relative  island  positions  must  be  correlated.  Thus  the  purpose  of  the 
multiple  scattering  calculations  will  be  to  correct  the  mean  field  theories 
for  the  contribution  of  the  correlation  to  the  effective  dielectric  constant 
of  the  medium. 

The  theory  is  based  on  the  concept  of  a  local  dielectric  constant  that 

fluctuates  around  an  average  value  eQ  that  is  independent  of  position.  In 

terms  of  operator  formalism  the  field  equation  can  be  written  in  terms  of  a 

2 

non-stochastic  operator  H  =  VxVx-K  and  a  stochastic  term 
2  0 

Ht  =  Kq  (6e/e0).  Where  KQ  is  the  propagation  constant  of  an  average  medium 
with  constant  eQ  and  6c  is  the  fluctuating  part  of  the  dielectric  constant. 
This  field  equation  is 

t  =  tQ  ♦  Go  Hx  t  (1) 

where  tQ  is  the  electric  field  that  results  from  solving  the  non-stochastic 
operator  and  £  is  the  resulting  scattered  field.  Since  we  are  in  the  presence 
of  a  random  medium  and  we  will  assume  that  an  effective  medium  dielectric 
constant  can  be  defined,  an  average  of  this  equation  can  be  performed  to 
obtain 

<2>  =  ?  ♦  G„  <H,  ?> 
o  o  . 

where  G  is  the  Green's  function  associated  with  2  .  In  order  to  be  able  to 
o  o 

obtain  useful  values  for  the  effective  dielectric  constant,  the  First 

6  7 

Smoothing  Approximation  ’  will  be  used.  Following  the  procedure  of  Karal  and 
s 

Keller  this  approximation  was  used  to  obtain  an  effective  medium  dielectric 
constant  c*°  this,  it  was  necessary  to  assume  a  plane  wave  solution 

for  the  macroscopic  field  with  an  effective  propagation  constant 

2  2  2 
K  3  U>  /c  £eff 


Previous  work  by  the  authors  indicate  that  the  pair  correlation  can  be 
expressed  as  an  exponential  relationship  with  a  correlation  length  L.  In  the 
spectral  range  of  interest  the  final  dispersion  relation  becomes 

i  ♦  2*1  ♦  ?£_-  [i  * 

2cq  3  a2  l-i2a  >-a  a3 
2  2  2  2 

where  a  =  eQw  L  /c  and  A  is  the  square  of  the  standard  deviation  of  the 
fluctuation. 

The  new  effective  dielectric  constant  includes  the  averaged  dielectric 
constant  plus  a  term  that  is  a  function  of  the  correlation  length  L.  This 
calculation  is  conceptually  based  on  the  assumption  that  there  is  an  averaged 
dielectric  constant  around  which  local  fluctuations  due  to  the  metallic  and 
dielectric  islands  occur,  known  as  the  Random  Continuum  Model.  Since  this  is 
precisely  the  model  that  is  described  by  the  Bruggeman  theory,  we  have  chosen 
as  the  averaged  dielectric  constant  eQ  the  one  obtained  from  that  theory. 

The  theoretical  predictions  of  both  the  Bruggeman  and  FSA  models  will  now 
be  compared  with  experimental  data  for  Ag/MgO,  80%  MgO,  and  Au/A1203,  84% 
A1203.  Figures  2  shows  the  imaginary  parts  of  the  index  of  refraction 
of  these  cermets.  For  the  calculation,  as  presented  in  the  Fig.  2a,  the 
dielectric  constant  of  the  metallic  islands  was  obtained  from  Orude's  model 
corrected  for  a  restricted  mean  free  patH  due  to  the  size  effects  and  due  to 
impurities  and  defects.  As  observed,  the  Bruggeman  model  predicts  higher 
values  for  the  imaginary  part  of  the  index  of  refraction  N2,  than  those 
obtained  experimentally,  whereas  the  FSA  model  fits  very  well  with 
experimental  data  in  almost  all  the  region  in  the  case  of  Ag/MgO  where  the 
radius  of  the  islands  is  of  the  order  of  80A  and  restricted  mean  free  path  of 
5A.  In  the  case  of  Au/A1203,  it  was  necessary  to  make  quantum  corrections  in 
order  to  obtain  the  best  fit,  because  the  islands  in  this  case  are  of  the 
order  of  30A.  In  this  second  case,  the  resonance  in  the  imaginary  part  of  N2 
is  shifted  towards  longer  wavelengths  with  respect  to  the  experimental 
results. 

Renormalized  FSA  Approximation  (RFSA) 

The  results  obtained  using  the  QSE  and  multiple  scattering  corrections  in 
the  FSA  approximation  to  the  B  model  shows  improvements  on  the  predicted  value 
for  the  imaginary  index  of  refraction,  but  these  corrections,  although 
affecting  the  magnitude  of  the  resonance,  do  not  broaden  or  shift  the 
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resonance  peak.  This  suggests  that  a  renormalized  approach  to  the  multiple 
scattering  calculation  may  be  the  correct  path  to  obtain  the  necessary 
corrections  to  the  position  and  width  of  the  resonance  peak.  A  renormalized 
procedure  can  be  devised  starting  from  the  field  Eq.  1  and  then  using  a 
T-matrix  calculation  for  single  spherical  scatterers.  Contrary  to  the  case  of 
the  FSA  approximation,  the  cermet  is  conceived  here  as  spherical  islands  of 
the  minority  constituent  embedded  in  a  matrix  of  the  dominant  constituent. 
Each  one  of  these  spherical  islands  is  then  taken  to  be  a  scattering  source 
from  which  the  overall  scattering  wave  for  the  random  medium  can  be  obtained. 
Since  this  model  of  cermet  is  consistent  with  the  MG  model,  the  whole 
procedure  will  be  constructed  on  a  zeroth  order  scattering  dielectric  constant 
eQ  that  is  obtained  from  a  MG  theory,  thus  assuming  a  non-correlated  medium  in 
which  the  spherical  grains  serve  as  scattering  sources.  Applying  now  the  First 
Smoothing  Approximation  yields  an  equation  that  is  formally  equivalent  to  the 
FSA  but  renormalized  and  written  in  the  following  form 

<!>  =  tQ  ♦  Gq'  <T  Gq'T>  <t>  ,  (2) 

where  G  '  is  the  renormalized  Green's  function  associated  with  <?  >. 
o  o 

The  second  term  in  the  r.h.s.  part  of  this  self-consistent  equation 
represents  scattering  processes  where  the  effective  total  field  <I>  is 
scattered  by  pairs  of  correlated  islands  a  and  B  located  respectively  at  r{ 
and  r..  In  order  to  apply  this  formalism  an  explicit  mathematical 

J 

representation  of  the  expression  in  the  second  term  of  the  r.h.s.  of  Eq.  (2) 

i  i 

must  be  obtained.  Peterson  and  Strom  developed  a  formalism  to  obtain  the 

T-matrix  for  an  array  of  n-scattering  objects  with  well-defined  spatial 

coordinates.  Their  formalism  is  an  extension  to  a  multiple  number  of 

„  12 

scatterers  of  the  one  developed  by  Waterman  for  a  single  scatterer.  In 
Waterman's  formalism  the  T-matrix  for  a  single  scatterer  is  obtained  expanding 
all  the  local  fields  in  terms  of  Spherical  Harmonics,  and  Bessel  or  Hankel 
functions.  The  term  describing  this  double  scattering  process  in  Peterson  and 
Strom's  formalism  is  given  by 

T  Gq'T  -  R< r j )  T*  o^-^)  TP  RCr^r.)  R(-f j)  ,  (3) 

where  T®,  T^  are  the  single  scatterer  T-matrices  of  islands  a  and  p, 

respectively,  expressed  with  their  center  in  the  origin  of  the  coordinate 
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Figure  2.  Comparison  between  experimental  data  10  and  calculations  for  Vj- 
(a)  Ag/WgO  (SOtMgO)  cermet  compared  with  B  (doted  line)  and  F$A 
(dashed  line)  models  with  1-125A,  and  a  restricted  mean  free 
path  «  IQA.  (b)  Au/Al^Oj  (845  aUOj)  cermet  compared  with  BO 
(doted  line),  80  (dashed  line)  and  FSAQ  (solid  line)  models  with 
l«80A,  a«30A  and  the  mean  free  path  »  1CA. 


Figure  3.  The  imaginary  part  Nj  of  the  refractive  index  is  drawn  for  an  Ag/*gO 
80%  MgO  cermet.  Solid  line  is  the  RF$A  calculation,  dashed  lines  that 
of  the  MG  model  and  dots  are  the  experimental,  (a)  n2  using  the  0 
model  with  t*  2.14x10  ^®$1n  both  MG  and  RFSA.  The  radius  of  the  is¬ 
lands  is  taken  as  80  A  and  l«  6  islands,  (b)  ?i2  using  QSE  with  a  re¬ 
stricted  mean  free  path  of  10  A  for  dielectric  constant  of  the  metallic 
grains  in  both  curves.  The  island  radius  is  125  A  and  l  same  as  (a). 
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system.  The  formalism  then  translates  the  particle  through  matrices  R  and  a 

to  their  correct  position  r,  and  r..  Expressions  for-T01  are  described  in 
11  12  1  J 

several  articles  ’  and  expressions  for  matrices  R  and  a  can  be  found  in 
Peterson  and  Strom's  article.  Using  Eq.  (3)  the  field  equation  can  now  be 
expressed  as  — 

<E(r)>  =  <E  (r)>  ♦  1/V2  2- 

cr ,  p*a 

x  Tp  R(r\-?.)  R(--?.)  <?>  d?.  dr. 

where  g  (r. ,F.)  is  the  two-point  correlation  function  and  V  is  the  volume  of 

*  J 

the  film.  We  again  consider  an  exponential  decay  type  correlation  function. 

To  simplify  the  calculation  we  assume  all  the  metallic  grains  to  be 

spherical  in  shape  and  of  equal  radius  which  allows  us  to  write  and  make 

them  diagonal  matrices.  The  field  <?>  and  <E  >  will  be  expanded  in  terms  of 

the  basis  functions  {Reij»  (K  r)}  that  are  solutions  to  the  Helmholtz  equation 

?2  o 

described  by  Waterman  KQ  is  the  propagation  constant  of  the 

electromagnetic  wave  in  the  zeroth  order  effective  medium.  Substituting  these 
expansions  for  the  fields  in  Eq.  (4)  it  is  possible  to  reduce  this  expression 
to  an  equivalent  set  of  scalar  equation 

Xn  *  d„  *  /  V  V  °nn'  d?j 

where  M  contains  all  the  correlation  effects, 
n 

The  coefficient  dn  representing  the  unperturbed  wave  can  be  written  in 
terms  of  the  effective  propagation  constant  K  ,  while  the  coefficient 
representing  the  total  average  field  <?>  can  also  be  expressed  in  terms  of  the 
effective  propagation  constant  K  of  the  correlated  medium  as: 


//  g(? 


g(r.,r*j)  R(rv)  f1  o(r.-r.) 


(4) 


dn  =  d®  exp[i?0*r]  ,  Xn  =  X®  exp[iK*f]  and  Bn  =  X®  exp[ i?*(r-r\)J 
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Following  the  procedure  described  by  Varadan,  et.  al.  we  can  use  the  scalar 
Helmholtz  operator  in  the  zeroth  order  approximation  to  simplify  the  above 
equation  and  reduce  the  volume  integral  to  a  surface  integral: 


where 
■nn'  = 


SV(K*-K02> 


*n  -  l  "n1  »„*■  'no'  •  (5> 

f  *° nn'(‘Ko?)  3exP<^,f?)/3r  ‘  exp(iK-f)  dann,  (-K0?)/8r}  ds. 

J|?|>a 
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Eq.  (5)  represents  a  system  of  homogeneous  coupled  equations  with  unknown 
coefficients.  Solving  the  associated  secular  equation  we  obtain  a  new 
dispersion  relationship  that  takes  into  account  pair  correlations  in  the 
scattering  medium.  In  terms  of  the  dielectric  constants  this  expression  can  be 
expressed  as 


e  =  r0  (1  +  2Fk2e0)/(l-Fk2F0)  ,  (6) 

where  e  is  the  new  effective  dielectric  constant  in  the  correlated  medium  of 
the  film  obtained  from  the  RFSA  model,  k  =  w/c  is  the  propagation  constant  of 
the  electromagnetic  wave  in  vacuum,  and  F  contains  all  the  effects  due  to 
correlations  between  islands  and  can  be  expressed  as 

F=  22/50  f2  exp(-2a/L)  (L2+2aL)  [(er-l)/(ep+2)]2 

where  f  is  the  volume  fraction  of  metal  and  e  is  the  ratio  of  the  dielectric 

r  ' 

constant  of  the  metallic  islands  to  the  external  medium. 

Results  and  Conclusions 

QSE  becomes  an  important  correction  to  the  effective  dielectric  constant 
of  the  metallic  grains  when  their  radius  is  less  than  100A.  In  using  the  QSE 
results  it  is  necessary  to  average  the  discrete  absortion  peaks  with  a 
log-normal  island  sizes  distribution,  before  calculating  the  optical  properties 
of  cermets.  The  next  step  in  improving  the  calculations  of  the  optical 

properties  is  to  include  two  point  correlation  between  neighboring  grains. 
In  order  to  simplify  the  multiple  scattering  correction,  the  First  Smoothing 
Approximation  was  used. 

In  performing  this  calculation  an  average  dielectric  constant  is  assumed 
that  is  conceptually  compatible  with  the  B  model  of  a  cermet.  Multiple 
scattering  corrections  in  the  FSA  approximation,  with  a  metallic  dielectric 
constant  corrected  for  QSE  in  the  cases  that  the  islands  are  smaller  than 
100A,  was  shown  to  improve  the  predicted  value  of  B  theory  but  was  unable  to 
significantly  shift  the  position  or  broaden  the  resonance  peak  of  the  complex 
index  of  refraction.  For  this  reason  it  was  considered  necessary  to  further 
improve  the  calculation  by  taking  into  account  renormalized  effects  in  a 
T-matrix  context.  The  result  of  this  calculation  shown  in  Eq.  6,  permits  the 
calculation  of  the  complex  effective  dielectric  constant  of  cermets  taking 
into  account  correlation  effects  that  include  renormalized  corrections. 


' .  * 


10 


We  have  applied  this  refined  calculation  to  Ag/MgO.  In  making  the 
calculations,  the  average  radius  of  the  islands  was  estimated  and  the  mean 
free  path  and  correlation  length  were  then  adjusted  to  obtain  the  best  fit 
with  experiment.  When  the  Orude  model  was  used  in  the  Ag/MgO  cermet  (See  Fig. 
3)  it  was  necessary  to  utilize  a  3A  mean  free  path  in  order  to  obtain  the  best 
fit  with  the  experimental  complex  imaginary  index  of  refraction.  This  was 
considered  to  be  too  short  a  mean  free  path  even  after  taking  into  account 
that  silver  tends  to  have  a  large  number  of  internal  imperfections  as 
mentioned  by  H.  Craighead10.  In  the  calculation  the  optimal  radius  for  the 
metallic  islands  was  80A,  a  value  considered  to  be  somewhat  small  according 
to  estimates  made  from  micrographs.  An  improvement  in  this  calculation  was 
obtained  when  the  average  radius  of  125A  for  the  islands  was  used  and  the  QSE 
corrections  were  made  (Fig.  3).  For  this  case,  the  shift  in  the  peak  was  good 
but  relatively  narrow  compared  to  the  experimental  results.  In  both  cases  the 
correlation  length  was  equal  to  6  islands,  a  value  that  is  compatible  with 
experimentally  measured  correlations  in  other  cermets. 

References 

1.  These  two  models  are  described  elsewhere.  The  original  papers  are:  3. 
C.  Maxwell -Garnett,  Phil.  Trans.  R.  Soc.  London  203  (1904)  385.  D.  A 
Bruggeman,  Ann.  Phys.  (Leipz)  24  (1935)  636. 

2.  R.  W.  Cohen,  G.  0.  Cody,  H.  0.  Coutts  and  B.  Abeles,  Phys.  Rev  B8  (1973) 
3689. 

3.  A.  Kawabata  and  R.  Kubo,  J.  Phys.  Soc.  J  pn.  21  (1966)  1765. 

4.  D.  M.  Wood  and  N.  W.  Ashcroft,  Phys.  Rev.  B25  (1982)  6255. 

5.  M.  Gomez,  G.  Rodriguez  and  L.  Fonseca,  Ferroelectrics  54  (1984)  223  and 
Ferroelectrics  54  (1984)  227. 

6.  R.  C.  Bourret,  Nuovo  Cimento  26  (1962)  1. 

7.  U.  Frisch  in  8harucha-Reid  (ed. )  Probabilistic  Methods  to  Applied 
Mathematics.  Academic  Press,  New  York,  1968,  p.  75. 

8.  F.  C.  Karal  and  J.  B.  Keller,  J.  Math.  Phys.  5  (1964)  537. 

9.  M.  Gomez,  L.  Fonseca  and  G.  Rodriguez,  Ferroelectrics  Lett. ,  2  (1984)  17. 

10.  H.  Craighead,  Ph.  D.  Thesis,  Cornell  University,  1980. 

11.  B.  Peterson  and  S.  Strom,  Phys.  Rev.  08  (1973)  3661. 

12.  P.C.  Waterman,  Phys.  Rev.  03  (1971)  825. 

13.  V.K.  Varadan,  V.N.  Bringi  and  V.V.  Varadan,  Phys.  Rev.  019  (1979)  2480. 


I 
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Multiple-scattering  theories  including  correlation  effects  to  obtain  the 
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Luis  Fonsecat,  Gerardo  Rodrigueztt,  Angel  Velazquez  and  Luis  Cruz.  University 
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Abstract 

Existing  theories  for  the  effective  optical  properties  of  a  ceramic-metal 
granular  medium  (cermet),  are  discussed  within  a  multiple  scattering  analysis 
and  used  as  the  starting  point  for  taking  into  account  correlation  effects. 
Analysis  of  the  nucleation  and  growth  process  in  cermets,  indicate  that  the 
correlation  between  the  islands  of  the  less  abundant  material  are  important  in 
the  study  of  optical  properties.  Multiple  scattering  corrections  are 
considered  by  solving  the  field  equation  in  the  First  Smoothing  Approximation. 
The  results  obtained  with  this  model  are  compared  with  other  theories  and 
against  experimental  data  of  real  cermets.  Quantum  size  effects  corrections 
to  the  dielectric  constant  of  the  metallic  grains  in  dielectric-1  ike  cermets 
and  their  effect  on  the  results  of  the  theory  are  also  analyzed.  Finally, 
renormalization  effects  introduced  through  a  T-matrix  formalism  are  discussed. 
Comparison  between  the  results  of  the  proposed  model  and  the  experimental  data 
for  different  cermets  shows  that  the  developed  theory  has  better  predictive 
value  than  earlier  models  and  good  agreement  is  obtained  with  experiment  even 
when  the  metal  concentration  is  large. 

tOn  leave  from  Escuela  de  Fisica,  Universidad  de  Costa  Rica. 
ttPresent  Address:  Cornell  University. 


I.  INTRODUCTION 


Cermets  are  inhomogeneous  materials  consisting  of  immiscibles  mixtures  of 
insulators  and  metals.  The  optical  properties  of  cermets  have  been 
extensively  researched  because  of  their  potential  application  to  the 

development  of  optically  selective  surfaces.1  Considerable  work  has  been 
performed  to  develop  theoretical  models  which  will  correctly  predict  the 
optical  properties  of  cermets  when  the  bulk  dielectric  constant  of  the 

constituents  is  known  experimentally.  It  has  been  found  that  the  optical 

properties  of  cermets  can  be  modified  by  changes  in  the  relative  concentration 
of  the  constituents  and  the  concomitant  change  in  microstructure. 2  Based  on 
constituent  concentration,  cermets  are  classified  into  three  general  types: 
a)  metallic  cermets  where  the  dominant  volume  fraction  is  metallic  and  as  a 
consequence  metallic  properties  dominate  the  optical  properties  of  the 

material;  b)  dielectric  cermets  where  the  dielectric  constituent  is  dominant, 
and  as  a  result,  metallic  inclusions  are  formed  in  a  dielectric  matrix; 
c)  intermediate  cermets  where  an  interconnected  metallic  network  develops,  in 

3 

this  region  percolation  effects  are  observed  .  At  opposite  extremes  when  the 
volume  fraction  of  the  metallic  or  dielectric  constituents  is  very  high 
(^•98%),  the  optical  properties  are  determined  by  the  major  constituent.  But 
when  the  concentration  of  the  dominant  constituent  is  less  than  98%  the 
material  develops  a  distinct  resonance  "which  is  not  characteristic  of  either 
one  of  the  constituents.  The  prediction  of  this  resonance  has  been  the  goal 
of  all  theories  developed  for  the  optical  properties  of  cermets  materials. 

Since  for  the  region  of  interest  in  the  electromagnetic  spectrum  the 

wavelength  of  the  radiation  is  large  compared  to  the  metallic  islands,  an 

electrostatic  approximation  can  be  made.  In  this  regime  a  simple  relation 

between  the  internal  (£.  )  and  the  external  (t  *)  electric  fields  for  a 

in  out 
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single  spherical  island  can  be  made  yielding  the  relationship 


2.  = 

in 


which  predicts  a  resonance  when  the  dielectric  constant  Ree(iu)  =  -2. 
Essentially  this  singularity  is  responsible  for  the  resonance  structure  in 
cermets. 

The  technological  usefulness  of  cermets  when  applied  to  optical  devices 
comes  from  their  optical  response  in  the  region  between  0.8pm  <  \  <  2.5pm. 
Since  in  the  spectral  region  of  interest  the  wavelength  is  larger  than  the 
microstructure,  an  effective  medium  dielectric  constant  can  be  used  to 
characterize  the  optical  response  of  the  materials.  Based  on  this  idea  several 
theoretical  models  have  been  proposed  to  obtain  the  effective  dielectric 
constant  of  cermets  in  terms  of  the  dielectric  constant  of  the  constituents, 
and  other  parameters  describing  the  microstructure,  such  as  island  radiuses, 
relative  concentrations  of  the  constituents,  and  the  shape  of  the  islands. 
Two  basic  and  widely  used  models  are  the  Maxwel 1 -Garnett  (MG)  and  the 
Bruggeman  (B)  models.5  The  MG  model  is  asymmetric,  since  it  considers  the 
cermet  as  composed  of  spheres  of  the  less  abundant  constituent  embedded  in  a 
matrix  of  the  most  abundant  material.  The  effective  dielectric  constant  is 
then  obtained  by  performing  a  volume  average  over  the  local  fields.  While  the 
B  model  is  symmetric  since  it  considers  the  two  constituent  materials  as 
spheres  embedded  in  an  effective  medium,  a  volume  average  is  also  performed  to 
obtain  a  self-consistent  expression  for  the  effective  dielectric  constant.  In 
both  theories  a  resonance  in  the  spectral  response  of  the  cermet  is  obtained, 
but  neither  agrees  with  the  experimentally  measured  complex  dielectric 
constant  (Fig.  1).  For  example,  comparison  between  the  experimentally 
obtained  imaginary  part  of  the  effective  index  of  refraction  and  the  MG  model 
reveals  that  the  theoretical  resonance  is  too  narrow  and  in  many  cases  blue 
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shifted  with  respect  to  experimental  results.  Several  papers,  like  those  of 
6 

Granqvist  show  that  the  magnitude  and  position  of  the  resonance  is  a  strong 
function  of  the  inclusion's  shape.  Corrections  to  the  dielectric  constant  of 
metallic  grains  for  size  effects  have  been  extensively  used  to  reduce  the 
magnitude  of  the  resonance,  but  they  fail  to  produce  enough  red  shift  and 
broadening  to  explain  experimental  data. 

The  objective  of  this  paper  will  be  to  show  that  corrections  due  to 

correlation  between  the  position  of  the  metallic  islands,  resulting  from  the 

formation  process  of  the  cermet,  can  move  and  broaden  resonance  obtained  by 

the  traditional  mean  field  theories  of  Maxwel 1 -Garnett  and  Bruggeman.  In 

order  to  perform  this  correction  a  multiple  scattering  theory  has  been 

developed  that  starts  from  an  effective  dielectric  constant  that  is  calculated 

assuming  no  correlation  between  the  metallic  islands  and  labeled  in  this  work 

as  the  zeroth  order  approximation.  Non  renormalized  and  renormalized  theories 

will  be  discussed  in  which  correlation  effects  are  taken  into  account.  The 

final  result  for  the  effective  dielectric  constant  will  be  expressed  in  terms 

of  the  zeroth  order  effective  dielectric  constant  e  (which  contains  no 

o 

correlation  effects  but  has  been  renormalized),  the  relative  concentration  of 
the  constituents  f,  the  dielectric  constant  of  the  less  abundant  constituent 
corrected  for  size  effects  em,  the  radius  of  the  inclusions  a,  and  a 
correlation  length  L. 

Existing  theories  including  those  that  take  into  account  the  shape  of  the 
islands  in  an  effort  to  improve  the  predicted  value  of  the  mean  field  theories 
will  be  analyzed  and  discussed.  Also  corrections  for  size  effects  on  the 
dielectric  constant  of  the  metallic  grains  will  be  studied  and  it  will  be 
shown  that  it  is  necessary  to  take  into  account  quantum  size  effects  for 
metallic  islands  of  good  metals  with  radiuses  less  than  100A. 
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A.  MEAN  FIELD  MODELS 

In  order  to  explain  the  characteristic  resonance  peak  of  cermets  and 
other  non-homogeneous  media  two  basic  mean  field  models  have  been  proposed  by 
J.  C.  Maxwel 1 -Garnett  and  Bruggeman.  Many  derivations  of  these  theories  have 
been  made  by  several  authors,7’8  and  it  has  been  shown  that  MG  theory  is  a 
zeroth  order  approximation  of  more  complex  multiple  scattering  theories.9’10 
The  authors  have  been  able  to  show  that  Bruggeman  can  also  be  reduced  to 
zeroth  order  approximation  in  a  multiple  scattering  formulation  of  the 
problem.  Both  theories  are  effective  medium  theories  within  the  electrostatic 
approximation  and  are  obtained  from  volume  averages.  Maxwel 1 -Garnett  theory 
is  an  asymmetric  model  considering  the  cermet  as  composed  of  islands  of  the 
less  abundant  constituent  embedded  in  a  matrix  of  the  most  abundant  material, 
while  the  Bruggeman  model  is  a  self-consistent  theory  and  assumes  that  both 
constituents  form  spheres  that  are  surrounded  by  the  effective  medium.  The 
response  of  islands  in  the  cermet  to  an  external  electric  field  produces  a 
resonance  in  the  complex  dielectric  constant.  In  general,  Maxwel 1 -Garnett 
produces  resonance  peaks  that  are  too  narrow  and  usually  shifted  towards  the 
blue  region  of  the  spectrum  with  respect  to  the  experimental  results,  while 
Bruggeman  tends  to  underemphasize  the  structure. 

Several  modifications  have  been  performed  in  an  attempt  to  obtain  better 
agreement  with  experimental  results.  In  particular,  considerable  effort  has 
been  invested  in  correcting  the  theory  for  effects  of  changes  in  the  shape  of 
the  islands  in  the  cermet.11'  Granqvist  6  has  taken  into  consideration  not 
only  shape  effects,  but  the  contribution  of  chains,  double  spheres,  fee 
clusters,  and  even  combinations  of  these  shapes  in  order  to  bring  mean  field 
theories  calculations  closer  to  experimental  results.  In  particular  he 
calculates  the  logarithm  of  the  transmission  for  metallic  type  cermet 
(Ag/Si02;  17%  SiOj)  using  the  MG  model,  under  the  assumption  that  the 
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dielectric  inclusions  are  ellipsoids  instead  of  spheres,  and  expressing  the 


effective  dielectric  constant  as: 


£MG  =  £ma 


1  +  2/3  for 
1  -  1/3  fa 


where  e  ,  is  the  bulk  dielectric  constant  of  the  matrix,  f  volume  fraction  of 
ma 

the  dielectric  and  a  is  proportional  to  the  polarizability.  Under  the 
assumption  that  ellipsoids  are  randomly  oriented,  the  value  for  a  is  given  by: 


..  -l  Z  — in  — 

3  i=l  e  +  D-  (e  •  ■  £  )  , 
ma  i  v  mi  ma' 

where  the  D-  are  the  triplets  of  the  depolarization  factors,  em-  is  the  bulk 
dielectric  constant  of  the  less  abundant  constituent,  and  the  1/3  term  is  a 
consequence  of  assuming  a  random  orientation  of  the  ellipsoids.  He  then  uses 
a  log-normal  distribution  for  the  ratio  of  the  major  semi-axis  a  to  the  minor 
semi-axis  c,  a  ratio  that  is  directly  connected  to  the  D's,  to  calculate  the 
effective  dielectric  constant  of  the  cermets.  Granqvist  centered  his 
log-normal  distribution  on  a  value  of  D  that  corresponds  to  the  ratio 
a/c  =  8.7,  with  this  ratio  he  weighted  the  distribution  towards  oblate 
ellipsoids.  Following  this  procedure  he  obtained  better  agreement  with 
experiment  than  when  all  the  inclusions  are  assumed  to  be  spheres. 

Figure  2  illustrates  the  dependence  of  the  imaginary  part  of  the 
dielectric  constant  e2  on  the  depolarization  factor  D  when  it  is  varied  from 
.2  <  D  <  1.0.  From  the  figure  it  can  be  seen  that  e-2  is  strongly  dependent 
on  the  depolarization  factor.  Several  distributions  were  used  by  us  to 
perform  the  average  with  the  depolarization  factors  assuming  values  between  .2 
and  .8,  values  that  were  considered  to  be  consistent  with  observed  micrographs. 


By  taking  a  log-normal  distribution  with  a  0  value  centered  at  an  average  D 
value  that  corresponds  to  a  ratio  a7c  =  1  (spheres),  it  was  shown  upon 
calculating  the  logarithm  of  the  transmission  that  the  result  corresponds  to 
the  one  obtained  if  only  spheres  would  have  been  used  for  the  metallic 
inclusions.  In  fact,  we  were  able  to  show  that  averaging  over  a  distribution 
of  shapes  was  not  sensitive  to  the  type  of  distribution  or  the  average  value 
of  the  a/c  ratio  used,  unless  a  highly  asymmetric  or  a  strongly  weighted 
towards  oblateness  distribution  was  utilized.  We  concluded  that  the  average 
polarizability  equals  the  values  obtained  for  spherical  grains  as  long  as 
distributions  consistent  with  observed  micrographs  are  used.  For  this  reason, 
it  was  decided  that  all  multiple  scattering  calculations  would  be  performed 
only  for  spherical  islands,  since  on  the  average,  due  to  the  random 
distribution  in  shape  and  orientation,  the  system  behaves  like  a  cermet  made 
up  of  spherical  islands. 

II.  SIZE  CORRECTION  TO  THE  DIELECTRIC  CONSTANT  OF  THE  METALLIC  INCLUSIONS 

A.  QUANTUM  SIZE  EFFECTS  (QSE) 

In  order  to  calculate  the  effective  dielectric  constant  of  cermets  it  is 
necessary,  in  general,  to  correct  the  bulk  dielectric  constant  of  the  minority 
constituent  for  size  effects.  Normally  the  insulator  has  a  dielectric 
constant  that  is  essentially  frequency  independent  in  the  spectral  range  of 
interest,  while  the  metallic  constituent"  has  a  strong  frequency  dependence  in 
the  same  region.  When  the  metal  is  the  minority  constituent,  its  bulk 
dielectric  constant  must  be  corrected  for  size  effects  when  the  bulk  mean-free 
path  becomes  comparable  to  the  size  of  the  grains.  Traditionally  Drude's 
model  together  with  experimentally  obtained  interband  contributions  is  used  to 
calculate  the  dielectric  constant  of  the  metal.  In  this  procedure,  the  mean 
free  path  of  the  electrons  in  the  metallic  island  is  corrected  for  boundary 
scattering  at  the  interface  between  the  island  and  the  external  medium.  Many 
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researchers  have  observed  that  in  order  to  obtain  a  reasonably  good  fit 
between  theory  and  experiment,  it  is  necessary  to  use  mean  free  paths  that  are 
considerably  smaller  than  the  ones  attributed  only  to  the  size  of  the  island. 
These  corrections  can  be  explained  as  the  combination  of  boundary  scattering 
at  the  metallic  grains  surface  and  scattering  inside  the  island  due  to 
internal  imperfections  and  faults.  Apparently,  the  contributions  from  these 
faults  and  imperfections  is  the  dominant  factor  in  determining  the  mean-free 
path  in  most  materials.  This  is  not  surprising  since  in  many  cases  the 
formation  process  is  based  on  the  coalescence  of  small  islands  to  form  larger 
ones,  a  process  that  should  produce  substantial  amounts  of  imperfections  in 
the  metallic  grains.  The  Drude  theory  with  these  corrections  is  satisfactory 
for  metallic  grains  of  good  metal  with  radiuses  of  more  than  100A.  However, 
when  the  size  of  the  metallic  particles  become  comparable  to  the  electron 
wavelength,  quantum  corrections  become  mandatory  if  a  realistic  dielectric 
constant  is  going  to  be  obtained. 

Until  recently  no  good  theoretical  models  existed  for  calculating 
dielectric  constant  of  metallic  particles  with  radius  of  the  order  of  100A  or 
less.13  To  solve  this  problem  several  researchers  have  used  the  Quantum  Box 
Method  within  the  Random  Phase  Approximation  to  calculate  the  dielectric 
constant  of  the  metallic  particles. 

14 

Kawabata  and  Kubo  have  corrected  the  dielectric  constant  for  size 
effects  assuming  that  the  spectrum  of  the  quantum  levels  form  essentially  a 
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continuum,  while  Wood  and  Ashcroft  carefully  reworked  the  Quantum  Size 
Method,  fully  taking  into  account  the  spectrum  resulting  from  the  box 
quantization  and  obtained  a  dielectric  constant  for  the  particles  that  can 
be  expressed  as 


Ree(x)  =  1  +  44  _a_  I  m2  (m  2  -  m2) 

n*  a0  m=l  ■' 


x  z  m'2[A2-(x2+r2)in-(-ll 

m'=l  9 


Ime(x)  =  44  _a_  r  I  m2  (m  2  -  m2) 

/tT  a0  x  m=l  c 


x  *  m'2[A2+x2+r2l[l-(-l)m+m' 


where 


0  =  A3  [(A2  -  x2  +  r2)2  +  4x2  r2];  A  =  m'2  -  m2 
x  =  2hw  ma2/h2/i2;  f  =  hm^/te^;  mc  =  INT(K^a/n) 


£ ^  is  the  Fermi  energy,  the  Fermi  wave  number,  a  is  the  size  of  the  grains, 

a0  the  Bohr  radius,  and  INT  the  integer  part  of  the  argument. 

Recently  this  formalism  has  been  applied  by  the  authors  of  this  article 

to  real  cermets  in  the  region  where  the  metallic  inclusions  in  the  cermets  are 
16 

smaller  than  100A  and  better  results  have  been  obtained  for  the  dielectric 
constant  than  when  Drude's  model  is  used. 

Figure  3  compares  the  imaginary  part  of  the  dielectric  constant  taking 
into  account  quantum  size  effects  following  Wood  and  Ashcroft  (QSE)  as 
expressed  in  Eq.  (1)  and  Kawabata  and  Kubo  model  (KK)  with  the  one  calculated 
using  Drude's  model.  In  the  calculation  of  the  Drude's  model  the  size  of  the 
silver  particles  are  taken  to  be  30A  and  no  interband  contributions  to  the 
dielectric  constant  have  been  included.  This  model  was  corrected  for  size 
effect  using  an  effective  relaxation  time  that  corrects  the  bulk  result  with 
boundary  scattering  through  the  relationship 
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1/T 


eff 


l/tb  +  Vf/a 


where  tb  is  the  bulk  relaxation  time,  is  the  Fermi  velocity  and  a  is  the 
radius  of  the  silver  grains.  The  multiple  peaks  appearing  in  the  QSE 
calculations  are  due  to  the  discrete  energy  spectrum  of  an  electron  in  a  box. 
These  peaks  tend  to  decrease  in  magnitude  as  one  moves  from  the  infrared  to 
the  visible  region  of  the  spectrum. 

Extensive  computer  calculations  were  made  to  compare  the  KK,  QSE,  and 
Drude's  model  for  good  metals  such  as  silver,  gold,  and  nickel  and  it  was 
found  that  in  all  cases  when  the  size  of  the  islands  exceeded  100A  the  three 
models  coincide  in  the  obtained  value  of  the  dielectric  constant  of  the 
metallic  inclusions,  provided  that  the  KK  model  is  corrected  for  a  geometric 
factor  of  tt/6. 17  For  that  reason  Orude's  formula,  which  is  the  simplest  one, 
will  be  preferred  for  metallic  islands  of  radiuses  larger  than  this  value. 

B.  APPLICATIONS  OF  THE  QSE  TO  REAL  CERMETS 

We  will  now  apply  the  QSE  result  to  the  calculation  of  the  optical 
properties  of  real  cermets  and  compare  the  results  with  those  obtained  using 
Drude's  model.  In  this  section  interband  contributions  to  the  dielectric 
constant  of  the  metal  grains  will  be  fully  taken  into  account. 18 
Experimentally  the  resonances  observed  in  the  dielectric  constant  of  an 
isolated  metal  grain  due  to  quantum  size  effect  are  not  observed  in  the 
effective  dielectric  constant  of  a  real  cermet.  The  disappearance  of  these 
peaks  is  due  to  the  fact  that  the  metallic  grains  in  a  real  cermet  do  not  have 
a  unique  size,  but  instead  have  a  distribution  of  sizes  as  can  be  verified  by 
micrographs  of  the  materials.  For  this  reason  the  dielectric  constant 
obtained  from  the  QSE  calculations  for  one  isolated  particle  were  averaged 
over  a  distribution  that  realistically  represented  the  variations  in  size  of 
actual  cermets. 


If  an  histogram  of  the  observed  radiuses  of  the  metallic  grains  is  made, 
a  distribution  of  the  form  illustrated  in  Fig.  4  is  obtained  showing  that 
the  distribution  is  skewed  towards  large  radiuses.  The  skewedness  is 
consistent  with  the  formation  process  of  the  inclusions,  in  which  large 
particles  are  formed  from  smaller  ones  by  coalescence.  Actually,  a  log-normal 
distribution  of  the  following  form  is  usually  used  to  fit  the  experimental 
results. 19 

f  (R)  =  (l/,/2n)  1/a  (exp  { - 1/2 [ ln( R/Rq )/ct] ^ }) 

where  Rq  is  the  average  radius  of  the  grains  and  a  is  a  measure  of  the  width 
of  the  distribution  and  is  given  by  a  =  £n  (cr^n ) .  Experimentally  the  value 
aSLn  1S  f°und  t0  be  between  1.1  and  1.5  for  most  cermets.20  Averaging  the 
dielectric  constant  of  one  metallic  grain  with  a  log-normal  distribution 
function,  an  effective  dielectric  constant  for  the  metallic  grains  in  the 
cermet  is  obtained.  This  average  QSE  result  is  then  used  to  calculate  the 
effective  medium  dielectric  constant  using  the  Bruggeman  mean  field,  and 
multiple  scattering  corrections  theories. 

In  Fig.  5  the  unaveraged  and  averaged  quantum  calculations  are  compared 
for  a  particle  of  average  radius  Rq  =  30  A.  As  seen  from  the  figure,  the 
effect  of  taking  the  average  with  the  Itfg-normal  distribution  is  to  smooth  out 
the  quantum  structure  of  a  single  particle  when  using  the  QSE  approximation. 
Since  experimentally  the  dielectric  function  of  a  cermet  is  found  to  be  a 
smooth  function  of  frequency,  this  is  the  correct  dielectric  constant  to  be 
used  for  the  metal  in  calculating  the  effective  dielectric  constant  of 
cermets.  For  radiuses  in  excess  of  100A  it  is  found,  as  expected,  that  the 
Drude  and  the  averaged  quantum  results  coincide. 


Before  applying  these  results  to  real  cermets  it  is  important  to  discuss 
the  relaxation  time  concept  as  used  in  the  Drude  and  the  QSE  methods.  Usually 
the  Drude  model  utilizes  an  effective  relaxation  time  that  corrects  the  bulk 
mean  free  path  of  the  electrons  for  collisions  with  the  grain  walls.  However, 
the  relaxation  time  is  not  well  understood  for  particles  of  radiuses  below 
100A.  According  to  Wood  and  Ashcroft,  surface  scattering  corrections  in  terms 
of  a  mean  free  path  concept  is  not  appropriate  when  the  presence  of  the 
surface  actually  determines  the  eigenstates  of  electrons  within  the  particle's 
volume.  The  mean  free  path  of  electrons  in  the  bulk  metal  for  the  case  of 
silver  is  of  the  order  of  100A.  While  for  smaller  metallic  particles  grown  by 
sputtering  processes  the  mean  free  path  is  reduced  due  to  imperfections  within 
the  metallic  grain  that  result  from  the  fabrication  process.  Small  metallic 
inclusions  in  cermets  are  known  to  grow  by  coalescence  of  smaller  inclusions 
which  diffuse  on  the  substrate,  and  as  a  result,  a  significant  number  of 
imperfections  are  trapped  within  the  metallic  grain  during  the  formation 
process.  To  take  into  account  the  effect  of  these  imperfections  we  will 
introduce  a  restricted  mean  free  path  smaller  than  the  bulk  mean  free  path 
when  applying  the  QSE  calculation  to  obtain  an  effective  complex  dielectric 
constant. 

We  now  compare  with  experimental  data  the  results  obtained  when  the  QSE 
dielectric  constant  is  introduced  into  the  calculation  of  the  effective 
dielectric  constant  of  the  cermet.  The  theoretical  calculations  will  be 
compared  with  the  measured  refractive  index  for  the  Au/A1203  and  Ag/MgO 
cermets  obtained  by  Craighead.21  This  experimental  data  was  selected  because 
it  is  obtained  from  cermets  with  metallic  grains  of  average  radiuses  less  than 
50A,  and  consequently,  the  QSE  method  should  contribute  significantly  to  the 
complex  index  of  refraction. 
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Figure  6(a)  compares  the  experimentally  measured  imaginary  part  of  the 
complex  index  of  refraction  for  Au/A1203  cermets  with  a  90%  volume  fraction  of 
the  dielectric  constituent,  with  the  Bruggeman's  theory  using  Drude's  model 
(BO)  and  Bruggeman's  theory  using  the  QSE  correction  (BQ).  The  average  radius 
for  the  calculations  was  taken  to  be  20A,  and  the  value  used  for  the  limited 
mean  free  path  was  also  20A.  As  can  be  observed  from  the  figure  the  BD  model 
exaggerates  the  magnitude  of  the  observed  resonance  in  the  imaginary  part  of 
the  index  of  refraction  N2.  The  BQ  model  greatly  improves  on  the  BD  for  N2  by 
reducing  the  magnitude  of  the  peak  and  approximating  much  better  the  actual 
values  of  N2  in  the  infrared  region  of  the  spectrum.  However,  in  this  case 
there  is  no  appreciable  shifting  of  the  position  of  the  peak.  Figure  6(b) 
shows  the  same  experimental  result  for  Au/A1203  cermet  but  now  compared  with 
Maxwel 1 -Garnett,  using  for  the  metal  the  Drude  dielectric  constant  (MGD)  and 
the  QSE  (MGQ)  values.  As  we  have  anticipated  MG  predicts  a  too  narrow  and 
generally  blue  shifted  peak.  In  this  case  QSE  also  reduces  the  magnitude  of 
the  peak,  and  a  measurable  red  shift  of  the  curve  and  a  broadening  effect  is 
observed. 

III.  MULTIPLE  SCATTERING  THEORIES 
A.  THE  FIRST  SMOOTHING  APPROXIMATION  (FSA) 

When  thin  films  of  cermets  are  grown  by  cosputtering  or  coevaporation 
methods,  nucleation  processes  are  responsible  for  the  formation  of  metallic  or 
dielectric  islands  in  the  matrix  of  the  corresponding  more  abundant 
constituent.  Since  the  nucleation  process  is  governed  by  diffusion  effects, 
relative  island  positions  must  be  correlated.  Several  authors  have  developed 
a  method  to  measure  this  correlation  and  demonstrated  its  existence.22  The 
Multiple  Scattering  Theories  developed  here  will  assume  that  the  islands  of 
the  minority  constituent  are  correlated  and  that  the  corresponding  correlation 
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length  is  sufficiently  large  to  make  important  contributions  to  the 
calculation.  Thus  the  purpose  of  the  multiple  scattering  calculations  will  be 
to  correct  the  mean  field  theories  for  the  contribution  of  the  correlation  to 
the  effective  dielectric  constant  of  the  medium.  It  will  be  shown  that  these 
effects  are  important  in  obtaining  a  model  with  predictive  value  for  the 
optical  response  of  cermets. 

The  theory  is  based  on  the  concept  of  a  local  dielectric  constant  that 
fluctuates  around  an  average  value  eQ  that  is  independent  of  position.  This 
position-dependent  dielectric  constant  can  be  written  as 

e(r,u»)  =  e0(w)  +  6e(r,uj)  , 

where  6e(r,u))  corresponds  to  the  spatially  fluctuating  term.  The  equation 
governing  the  electric  field  of  the  electromagnetic  wave  can  then  be  written 
in  terms  of  this  spatial ly-dependent  dielectric  constant  as 

V  x  V  x  §  -  (uj/c)2pet  -  pVpxVx£  =  0 

Rewriting  this  equation  explicitly  in  terms  of  the  fluctuating  part  of  the 
dielectric  constant  it  can  be  expressed  as 

V  X  V  X  i  -  Ko2!  =  Kg  (6e/e0)  t 

where  K0  is  the  average  propagation  in  a  medium  with  average  dielectric 
constant  eQ.  In  terms  of  operator  formalism  the  field  equation  can  be  defined 

as  a  non-stochastic  operator 

2 

HsVxVx-K 

o 


and  a  stochastic  term 


2 

H  =  Ko  (5e/e0) 

The  field  equation  then  becomes 

N  ?o  +  Go  Ht  ?  ,  (2) 

where  i:  is  the  electric  field  that  results  from  solving  the  non-stochastic 
operator.  Since  we  are  in  the  presence  of  a  random  medium  and  we  will  assume 
that  an  effective  medium  dielectric  constant  can  be  defined,  an  average  of 
this  equation  can  be  performed  to  obtain 

<(:>  =  to  +  Go  <Hi  (■> 

From  this  expression  an  infinite  multiple  scattering  expansion  in  terms  of  Hj 
can  be  performed  with  the  condition  that  the  fluctuating  term  is  centered 
around  the  averaged  dielectric  constant  and  therefore,  cH^  =  0. 

Diagrammatical ly  this  expansion  can  be  represented  as  follows 


where  --  represents  Green's-  function  associated  with  j-  r~~~l  represents 
Green's  function  associated  with  <£>  and  0  represents  Ht.  The  dashed  line 
indicates  correlation  between  local  scattering  points.  The  above  equation 
although  exact,  cannot  be  solved.  In  order  to  be  able  to  obtain  useful  values 
for  the  effective  dielectric  constant  the  First  Smoothing  Approximation,23’24 
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which  only  takes  into  consideration  pair  correlation  functions,  will  be  used. 
This  approximation  yields  the  equation 

f:  =  £0  +  Go  <Ht  Go  Ht>  <(:>  , 

which  upon  interation  can  be  expressed  diagrammatical ly  as 

1 — 1  —  —  +  — 0 — b —  +  — d — b — 6 — 0 —  +  •  •  •  • 

This  diagram  contains  the  pair  correlation  functions 

<6  eCrJd  e(r2)> 

represented  by  the  dashed  lines.  Following  the  procedure  of  Karal  and  Keller25 
this  approximation  was  used  to  obtain  an  effective  medium  dielectric  constant 
e  ff.  To  do  this,  it  was  necessary  to  assume  a  plane  wave  solution  for  the 
macroscopic  field  with  an  effective  propagation  constant 

2  22 

K  =  (u)  /c  )  £g ^ 

Previous  work  by  the  authors26  indicate  that  the  pair  correlation  can  be 
expressed  as  an  exponential  relationship  with  a  correlation  length  L  in  the 
following  form 

<6£(r1)6e(r2)>  =  &  e  r/^L 

r  =p2  “  n|  and  A  was  shown  to  be  well  represented  by  the  square  of  the 
standard  deviation  of  the  fluctuations.  In  the  spectral  range  of  interest  the 
final  dispersion  relation  becomes 
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!*,f  =  £o  *  270 


where 


1  +  ill 1  + 


l-i2a 


-  ♦  is!  cot -lllsl± 


a2  =  e0w2L2/c2. 

The  new  effective  dielectric  constant  includes  the  averaged  dielectric 
constant  plus  a  term  that  is  a  function  of  the  correlation  length  L.  This 
calculation  is  conceptually  based  on  the  assumption  that  there  is  an  averaged 
dielectric  constant  around  which  local  fluctuations  due  to  the  metallic  and 
dielectric  islands  occur  and  is  known  as  Random  Continuum  model.  Since  this 
is  precisely  the  model  that  is  described  by  Bruggeman1 s  theory,  we  have  chosen 
as  the  averaged  dielectric  constant  eQ  the  one  obtained  from  that  theory. 
This  is  perfectly  consistent  with  multiple  scattering  theory  since  it  can  be 
shown  that  the  Bruggeman  model  is  a  zeroth  order  approximation  of  multiple 
scattering  theories.  The  proof  consists  in  assuming  that  both  the  metallic 
and  the  dielectric  constituents  form  spherical  scattering  sources  in  an 
average  medium,  and  then,  in  the  absence  correlation,  the  average  of  the 
T-matrix  <T>  can  be  written  as  the  volume  fraction  average  of  the  single 
particle  metal  T-matrix  Tt  and  the  dielectric  particle  T-matrix  T2  as: 
<T>  =  f  T1  +  (l-f)T2  where  f  is  the  volume  fraction.  Then  taking  the  long 
wavelength  limit  and  assuming  the  self-consistent  condition,  <T>  =  0,  the 
Bruggeman  equation  is  immediately  obtained. 

The  theoretical  predictions  of  both  the  Bruggeman  and  FSA  models  will 
now  be  compared  with  experimental  data  for  Ag/MgO;  80%  MgO,  and  Au/A1203;  84% 
A1203.  Figures  7(a)  and  7(b)  show  the  imaginary  parts  of  the  index  of 
refraction  of  these  cermets.  For  the  calculation,  as  presented  in  the 
Fig.  7(a),  the  dielectric  constant  of  the  metallic  islands  was  obtained  from 
Orude's  model  corrected  for  a  restricted  mean  free  path  due  to  size  effects 
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and  due  to  impurities  and  defects.  As  observed,  the  Bruggeman  model  predicts 
higher  values  for  the  imaginary  part  of  the  index  of  refraction  N2,  than  those 
obtained  experimentally,  whereas  the  FSA  model  fits  very  well  with 
experimental  data  in  almost  all  the  region  in  the  case  of  Ag/MgO  where  the 
radius  of  the  islands  is  of  the  order  of  80A  and  restricted  mean  free  path  of 
5A  is  used.  In  the  case  of  Au/A1203,  it  was  necessary  to  make  quantum 
corrections  in  order  to  obtain  the  best  fit,  because  the  islands  in  this  case 
are  of  the  order  of  30A.  In  this  second  case,  the  resonance  in  the  imaginary 
part  of  N2  is  shifted  towards  longer  wavelengths  with  respect  to  the 
experimental  results. 

B.  RENORMALIZED  FSA  APPROXIMATION  (RFSA) 

The  results  obtained  using  the  QSE  and  multiple  scattering  corrections  in 
the  FSA  approximation  to  the  B  model  shows  improvements  on  the  predicted  value 
for  the  imaginary  index  of  refraction,  but  these  corrections,  although  they 
affect  the  magnitude  of  the  resonance,  do  not  broaden  or  shift  the  resonance 
peak.  In  general,  the  B  and  MG  model  do  not  predict  correctly  the  position  or 
width  of  the  absorption  peak.  This  suggests  that  a  renormalized  approach  to 
the  multiple  scattering  calculation  may  be  the  correct  path  to  obtain  the 
necessary  corrections  to  the  position  and  width  of  the  resonance  peak.  A 
renormalized  procedure  can  be  devised24  starting  from  the  field  Eq.(3)  and 
then  using  a  T-matrix  calculation  for  sfngle  spherical  scatterers.  Contrary 
to  the  case  of  the  FSA  approximation,  the  cermet  is  conceived  here  as 
spherical  islands  of  the  minority  constituent  embedded  in  a  matrix  of  the 
dominant  constituent.  Each  one  of  these  spherical  islands  is  then  taken  to  be 
a  scattering  source  from  which  the  overall  scattering  wave  for  the  random 
medium  can  be  obtained.  Since  this  model  of  cermet  is  consistent  with  the  MG 
model,  the  whole  procedure  will  be  constructed  on  a  zeroth  order  scattering 
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dielectric  constant  that  is  obtained  from  a  MG  theory,  thus  assuming  a 
non-correlated  medium  in  which  the  spherical  grains  serve  as  scattering 
sources.  Eq.  (2)  can  be  written  in  diagrammatic  form  as 

a  a  0  a  fl  a 

G  =  —  +  — 0 —  +  — 0 — 0 —  +  —0—0—0—  +  •••  ,  (4) 

where  the  superscripts  on  each  point  refer  to  a  specific  scatterer  that 
produces  a  local  fluctuation  in  the  background  medium.  Defining  the  single 
scattering  T-matrix  as 

Ta  K  *  H>  «  t 

or 

a  a  a  a  a  a 
T  =  0  +  0-0  +  0-0-0  +**.=• 

Of 

where  the  symbol  0  represents  the  new  renormalized  potential  for  single 
scattering  spherical  sources.  In  terms  of  this  new  renormalized  potential 
Eq.  (4)  becomes 

G  =  —  +  —8—  +  --8—8--  +  --8—0—8—  +  ••• 
or  a  0  a  D  Y 

Naturally,  in  this  equation  the  subscripts  cannot  be  identical  in  two 
consecutive  points  since  these  terms  have  already  been  taken  into  account  in 
the  single  particle  scattering  T-matrix.  If  the  statistical  average  is 
performed  on  this  Green's  function,  the  following  statistically  averaged 
Green's  function  is  obtained 
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_  ^*V 

l  I  =  —  +  —ft—  +  +  —4— e--  +  ••• 

a  a  p  a  p 

this  expression  contains  both  correlated  and  uncorrelated  terms.  A 
renormalized  propagator  that  only  takes  into  account  uncorrelated  terms  is 
then  defined  as 

G  '  =  -n.  =  —  +  —8—  +  --a—*--  +  —8—8—8—  +  ••• 
a  op  a  p  a 

where  G  '  is  the  renormalized  uncorrelated  Green's  function  identified  with 
o 

the  propagation  of  an  electromagnetic  wave  in  an  uncorrelated  medium.  In 
terms  of  the  renormalized  Green's  function  and  scattering  potentials  the 
general  correlated  Green's  function  can  be  expressed  as 

*  + 

a  p  a  P  a 

This  expression  is  formally  equivalent  to  the  Random  Continuum  final 
Eq.  (3)  that  was  used  in  the  FSA  approximation  with  the  difference  that  now 
the  scattering  sources  and  the  intermediate  propagation  Green's  functions 
have  been  renormalized  taking  into  account  all  uncorrelated  terms.  A  simple 
correspondence  can  be  made  between  the  previous  and  the  present  equation 
through  the  following  symbolic  equivalences 

-■>*'«  that  is  G„  -*•  G  ' 
o  o 

0  ft  that  is  Ht  -*  T 

As  previously  stated,  in  the  absence  of  correlation  effects  or  equivalently 
in  an  uncorrelated  medium,  a  zeroth  dielectric  constant  is  obtained  that 


corresponds  to  the  MG  model.10  Applying  now  the  First  Smoothing  Approximation 
to  this  diagram,  which  is  equivalent  to  neglecting  all  correlations  higher 
than  the  two  point  correlation  and  also  neglecting  unlinked  and  nested  two 
point  correlations,  this  yields  an  equation  that  is  formally  equivalent  to  the 
FSA  but  renormalized  and  written  in  the  following  form 


<&  =  t  +  G  '  <T  G  'T>  <t> 

OO  0 


diagrammatical ly  this  can  be  written  as 


i>  =  <tQ>  +  <&— — ®  <t> 

a ,  0*or  a 


where  <tQ>  is  the  statistically  averaged  field  propagating  in  the  uncorrelated 
medium  and  <£>  the  effective  total  field.  The  second  term  in  the  r.  h.s.  of 
this  self-consistent  equation  represents  scattering  processes  where  the 
effective  total  field  <£>  is  scattered  by  pairs  of  correlated  islands  a  and  B 
located  respectively  at  r.  and  r..  In  order  to  apply  this  formalism  an 

J 

explicit  mathematical  representation  of  the  diagram  in  the  second  term  of  the 
r.h.s.  of  Eq.  (5)  must  be  obtained.  Peterson  and  Strom27  developed  a 
formalism  to  obtain  the  T-matrix  for  ah  array  of  n-scattering  objects  with 
well-defined  spatial  coordinates.  Their  formalism  is  an  extension  to  a 
multiple  number  of  scatterers'  of  the  one  developed  by  Waterman28  for  a  single 
scatterer.  In  Waterman's  formalism  the  T-matrix  for  a  single  scatterer  is 
obtained  expanding  all  the  local  fields  in  terms  of  Spherical  Harmonics,  and 
Bessel  or  Hankel  functions.  The  term  describing  this  double  scattering  process 
in  Peterson  and  Strom's  formalism  is  given  by: 
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"*■  R(^j)  RCr^-r^)  R(-rj) 


where  1“,  are  the  single  scatterer  T-matrices  of  islands  a  and  p, 
respectively,  expressed  with  their  center  in  the  origin  of  the  coordinate 
system.  The  formalism  then  translates  the  particle  through  matrices  R  and  a 
to  their  correct  position  r.  and  r..  Expressions  for  Ta  are  described  in 

J 

several  articles*7 ,28and  expressions  for  matrices  R  and  a  can  be  found  in 
Peterson  and  Strom's  *7  article.  Using  Eq.  (2)  the  field  equation  can  now  be 
expressed  as: 


d:(r)>  =  <£Q(r)>  +  1/V  ^  f  f  gO^.fj)  R(fj)  T° 

a,P*ofJ  J 


x  o(r.-rj)  Tr  RCr^.-r.)  R(-r ^ )  <E>  dr.  dr^ 

where  g  (r.,rj)  is  the  two-point  correlation  function  and  V  is  the  volume  of 
the  film.  Considering  the  spheres  as  impenetrable  and  the  medium  as  isotropic, 
an  exponential  metal-metal  correlation  function  can  be  defined  as: 

g^.rj)  =  0(r-2a)  exp  (-  |r.  -  r^J  /  L) 

where  0  is  the  Heaviside  function,  L  the  correlation  length,  and  the 
correlation  function  measures  the  excess  of  the  metal  density  with  respect  to 
its  average,  for  this  reason  the  function  goes  to  zero  at  large  distances.  To 
simplify  the  calculation  we  assume  all  the  metallic  grains  to  be  spherical  in 
shape  and  of  equal  radius  which  allows  us  to  write  Ta=T^  and  make  them 
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diagonal  matrices.  Assuming  a  macroscopic  volume  we  can  change  (1/V  )Z  by  6  , 
where  5  is  the  density  of  islands  in  the  cermet.  With  this  approximation  the 
field  equation  can  be  written  as 

<E(r)>  =  <EQ(r)>  +  6*  j*  drj  R(r.)Jdl  exp[-l/L]  T 

l>2a 

x  o(t)  T  R(-t)  R(-r.)  <t> 

where  the  substitution  t  =  r.  -  r.  has  already  been  performed.  The  field  <?> 
and  <IQ>  will  be  expanded  in  terms  of  the  basis  functions  { Re*^n(KQr) }  that  are 
solutions  to  the  Helmholtz  equation  described  by  Waterman.28  The  $n  are  known 
as  the  elementary  fields  and  KQ  is  the  propagation  constant  of  the 
electromagnetic  wave  in  the  zeroth  order  effective  medium.  In  terms  of  this 
base  the  fields  can  be  expressed  as: 

<t0>  =  Zdn  r4„(k0;) 

and 

<£(*•)>  =  IXn  Re$n(KQr) 


While  the  total  effective  field  <?>  in  the  second  term  of  the  r.h.s.  is 
expanded  in  a  base  centered  at  one  of  the  correlated  islands: 

<h  * 

Substituting  these  expansions  for  the  fields  in  Eq.  (7)  and  performing  the 
integration  over  the  relative  coordinate  between  correlated  islands  we  obtain 


(8) 


2*„  r4b<k0?>  =  2dn  R«Hn(K0?) 

*  Mn  {  Bn  d?j 

Matrix  elements  Mn  contain  all  the  correlations  effects  and  can  be 
written  in  terms  of  the  elements  of  the  single  scatterer  T-matrix.  Since  the 
wave  length  of  incident  radiation  relevant  to  our  problem  is  \  >  .3  pm  and  the 
radiuses  of  islands  for  typical  cermets  are  of  the  order  of  .01  pm,  the  long 
wave  limit  can  be  utilized,  thus  reducing  the  expansion  of  the  electromagnetic 
field  to  the  magnetic  and  electric  dipolar  terms.  In  this  limit  only  the 
first  six  terms  of  these  matrices  are  necessary: 

Mi  =  M3  =  M5  =  Ai  Ti  +  A2  Tj  T2 

M2  =  M4  *  M6  »  Ai  T22  +  A2  Ti  T2  ; 

where, 

Ai  s  -  22/ii/5K0  (L2+2al)  exp[-2a/L] 

A2  s  5/11  Ax 

Ti  s  i  2/15  (K0a)5  (l-er) 

T2  s  i  2/3  (K0a)3  (er-l)/(e/2) 

e  is  the  ratio  of  the  dielectric  constants  of  the  metallic  islands  and  the 
r 

external  medium,  and  it  represents  the  source  of  the  multiple  scattering 
processes.  To  perform  the  integral  over  the  r.  coordinate  the  properties  of 

J 

the  translation  matrices  o29  must  be  used  to  displace  the  $n  from  position  r. 
to  the  origin: 


then  Eq.  (8)  reduces  to 


?Xn  *  1%  RA<V> 

*525Mn/Bn  WKoVRein'(K o?)  dfj  ■  <9> 

All  terms  in  Eq.  (9)  are  now  written  in  the  same  orthogonal  base  and, 
therefore,  can  be  reduced  to  an  equivalent  set  of  scalar  equations: 

X  =  d  +  622  /  M  ,  B  ,  a  ,  dr. 
n  n  ~  /  n  n  nn  j 

The  limits  of  integration  in  the  above  expression  must  be  restricted  to  values 
|r j-r| >  a  in  order  to  avoid  having  the  scattered  field  inside  the  scattering 
island  located  at  r..  Since  in  the  long  wave  limit  the  electromagnetic  wave 

J 

tends  to  average  the  microstructure  at  the  scale  of  island  sizes  it  can  be 
assumed  that  the  wave  propagates  in  the  effective  medium  as  a  plane  wave.  The 
coefficient  dn  representing  the  unperturbed  wave  can  be  written  in  terms  of 
the  effective  propagation  constant  KQ  as: 

dn  =  d®  exp[iK0-r] 

While  the  coefficient  representing  the  total  average  field  <E>  can  also  be 
expressed  in  terms  of  the  effective  propagation  constant  K  of  the  correlated 
medium  as: 


0 


8n  =  Xn  expCiK-Cr—r^)] 


In  terms  of  these  expressions  the  field  equation  now  becomes: 


X*  exp(i£*r)  =  d®  exp(itfQ-r) 


*  Mn'  Xn'  [  «„„.«0V  dfj 

-*1+  -V  | 


r'rj|  ’  3 


Following  the  procedure  described  by  Varadan  et  al.  the  scalar  Helmholtz 
operator  can  be  used  in  the  zeroth  order  approximation  to  simplify  the  above 
equation  and  reduce  the  volume  integral  to  a  surface  integral: 

Xn=  ?H„.  X*  I„„,  ,  (10) 


where 


2  2  2 
‘nn'  '  6  ‘Ko  : 


■o  S)/j0m.l(V 
J  r  >  a 


)  3exp(iK-r)/8r 


exp(ift*r)  9ann' ("K0r)/8r]ds. 


Eq.  (10)  represents  a  system  of  homogeneous  coupled  equations  with  unknown 
coefficients.  Solving  the  associated  secular  equation  a  new  dispersion 
relationship  was  obtained  that  took  into  account  pair  correlations  in  the 
scattering  medium: 


«2  -  Kq2  =  (-i7t62M6/K0)[2  +  ( K/KQ ) 2  ] 


•*.  •*.  •*.  -• .  -* 


In  terms  of  the  dielectric  constants  this  expression  can  be  expressed  as: 

e  =  eQ  (1  +  2Fk2eo)/(l-Fk2£o)  ,  (11) 

where  e  is  the  new  effective  dielectric  constant  in  the  correlated  medium  of 
the  film  obtained  from  the  RFSA  model,  k  =  uj/c  is  the  propagation  constant  of 
the  electromagnetic  wave  in  vacuum,  and  F  contains  all  the  effects  due  to 
correlations  between  islands  and  its  explicit  form  is, 

F=  22/50  f2  exp(“2a/L)  (L2+2aL)  [(er~l)/(er+2)]2 
IV.  RESULTS  AND  CONCLUSIONS 

The  classical  approach  to  the  theory  of  the  optical  properties  of  cermets 
has  traditionally  used  Orude's  model  corrected  for  interband  transitions  to 
obtain  the  complex  dielectric  constant  of  the  metallic  constituent  with  an 
effective  relaxation  time  that  takes  into  account  surface  scattering  but  fails 
to  consider  quantum  effects.  From  the  present  work  it  is  evident  that  this 
approach  is  valid  when  the  radiuses  of  the  metallic  inclusions  are  larger  than 
100A  when  the  metals  used  are  good  metals.  We  have  shown  that  when  the 
radiuses  of  the  metallic  grains  are  smaller  than  100A,  QSE  becomes  an 
important  correction  to  the  effective  -dielectric  constant  of  the  metallic 
grains  when  calculating  the  optical  properties  of  cermets.  When  using  the  QSE 
results  it  is  necessary  to  average  the  discrete  absortion  peaks  with  a 
log-normal  distribution  of  island  sizes.  This  average  is  necessary,  since 
experimentally  no  structure  is  observed  in  obtained  optical  constants  of 
cermets.  In  general,  the  ability  to  reproduce  the  experimentally  observed 
optical  properties  of  cermets  is  considerably  improved  when  the  dielectric 


constant  of  the  metallic  constituents  for  islands  less  than  100A  is  corrected 
for  quantum  size  effects. 

The  next  step  to  improve  the  calculation  of  the  optical  properties  of 
cermets  is  to  include  two  point  correlation  between  neighboring  grains.  Under 
this  condition  a  multiple  scattering  theory  must  be  developed  that  will  take 
into  account  interference  effects  between  the  scattering  produced  by 
neighboring  metallic  grains  that  are  correlated.  In  order  to  simplify  the 
multiple  scattering  correction,  the  First  Smoothing  Approximation  was  used, 
thus  limiting  the  calculation  to  pair  correlations  through  a  self-consistent 
equation.  In  performing  this  calculation  an  average  dielectric  constant  is 
assumed  that  is  conceptually  compatible  with  the  B  model  of  a  cermet. 
Multiple  scattering  corrections  in  the  FSA  approximation,  with  a  metallic 
dielectric  constant  corrected  for  QSE  in  the  cases  that  the  islands  are 
smaller  than  100A,  was  shown  to  improve  the  predicted  value  of  B  theory  but 
was  unable  to  significantly  shift  the  position  or  broaden  the  resonance  peak 
of  the  complex  index  of  refraction.  For  this  reason,  it  was  considered 
necessary  to  further  improve  the  calculation  by  taking  into  account 
renormalized  effects  in  a  T-matrix  context.  The  result  of  this  calculation 
shown  in  Eq.  (11),  permits  the  calculation  of  the  complex  effective  dielectric 
constant  of  cermets  taking  into  account  correlation  effects  that  include 
renormalized  corrections.  This  equation  is  expressed  in  terms  of  the 
following  basic  physical  parameters  obtainable  from  experimental  data:  the 
average  radius  of  the  metallic  island  a;  the  correlation  length  L;  volume 
fraction  of  the  metallic  constituents  f;  the  bulk  dielectric  constant  of  the 
metallic  grains  and  the  bulk  dielectric  constant  of  the  dielectric 
material.  As  the  metallic  concentration  is  reduced  the  correlation  length 
tends  to  vanish  and  we  recover  from  Eq.  (11)  the  zeroth  order  approximation  e0, 
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the  effective  dielectric  constant  of  the  uncorrelated  medium.  Since  in  the 
renormalized  calculation,  islands  of  the  minority  constituent  are  assumed  to 
be  the  scattering  sources  embedded  in  a  medium  of  the  majority  constituent, 
the  proper  starting  zeroth  order  approximation  dielectric  constant  is  the  MG 
model . 

As  previously  stated  the  bulk  metallic  dielectric  constant  must  be 
corrected  for  size  effects  for  islands  of  100A  or  less.  We  have  applied  this 
refined  calculation  to  Ag/MgO  and  Ni /A1 2O3  cermets  where  experimental  data  is 
available  from  Craighead's  work.21  The  percentage  of  metallic  constituent  is 
20%  in  Ag/MgO  and  between  39  and  46%  in  N i /A1 2O3 .  These  concentrations  are 
considered  to  be  high  and  are  normally  not  predicted  well  by  mean  field 
theories.  Since  as  metal  concentration  increases  correlation  effects  become 
more  important,  it  is  expected  that  multiple  scattering  effects  should  be  more 
significant  in  these  high  concentration  regimes.  In  making  the  calculations, 
the  average  radius  of  the  islands  was  estimated  and  the  mean  free  path  and 
correlation  length  were  then  adjusted  to  obtain  the  best  fit  with  experiment. 
When  the  Drude  model  was  used  in  the  Ag/MgO  cermet  (Fig.  8(a))  it  was 
necessary  to  utilize  a  3A  mean  free  path  in  order  to  obtain  the  best  fit  with 
the  experimental  complex  imaginary  index  of  refraction.  This  was  considered 
to  be  too  short  a  mean  free  path  even  after  taking  into  account  that  silver 
tends  to  have  a  large  number  of  internal  imperfections  as  mentioned  by  H. 
Craighead. 21  In  the  calculation,  the  optimal  radius  for  the  metallic  islands 
was  80A,  a  value  considered  to  be  somewhat  small  according  to  estimates  made 
from  micrographs.  An  improvement  in  this  calculation  was  obtained  when  the 
average  radius  of  125A  for  the  islands  was  used  and  the  QSE  corrections  were 
made  (Fig.  8b).  For  this  case,  the  shift  in  the  peak  was  in  good  agreement 
but  relatively  narrow  compared  to  the  experimental  results.  In  both  cases  the 
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correlation  length  was  equal  to  6  islands,  a  value  that  is  compatible  with 
experimentally  measured  correlations  in  other  cermets. 

The  main  conclusion  that  can  be  obtained  from  this  calculation  is  that 
our  RFSA  theory  make  the  best  correction  to  the  MG  mean  field  theories.  It  is 
possible  to  play  around  with  the  three  adjustable  parameters  to  obtain  better 
fits,  but  what  is  significant,  is  that  all  values  are  compatible  with 
experimentally  known  facts. 

Figure  9  is  a  similar  comparison  for  Ni/Alj,03  with  metallic 
concentrations  of  39  and  46%.  Very  good  agreement  for  the  position  of  the 
resonance  peak  and  overall  width  was  achieved  with  reasonable  values  for  the 
mean  free  path  of  10A,  and  with  radiuses  consistent  with  the  micrographs. 21 
The  fit  between  theory  and  experiment  obtained  in  this  figure  is  particularly 
significant,  since  the  metallic  concentrations  used  are  close  to,  if  not 
within  the  percolation  regime,  a  region  where  the  mean  field  models  totally 
fail  to  reproduce  experimental  results.  Probably  this  regime  is  close  to  the 
limit  of  concentrations  for  which  the  RFSA  approximation  will  be  valid  since 
we  do  not  expect  our  calculations  to  be  applicable  in  the  percolation  region. 

The  results  of  this  and  other  recent  paperZB  leads  us  to  conclude  that 
correlation  effects  and  multiple  scattering  corrections  make  important 
contributions  to  the  optical  properties  of  cermet  materials  and  need  to  be 
included  in  theories  with  predictive  value.  Equation  (11)  is  restricted  to 
spherical  shapes  in  the  long  wave  limit,  but  the  formalism  as  presented  here, 
can  be  extended  to  include  other  island  shapes  and  to  permit  corrections  for 
higher  multipolar  contributions  to  the  field  equation.  These  multipolar 
corrections  may  become  important  when  the  metal  concentration  increases  to  the 
point  that  island  proximity  is  no  longer  compatible  with  the  dipolar 
approximation.  The  proposed  dispersion  relations  depend  strongly  on  the  model 


used  for  the  zeroth  order  dielectric  constant  and  the  values  used  for  the 
dielectric  constant  of  the  metallic  grains.  Other  zeroth  order  models  for  the 
effective  medium  dielectric  constant  are  being  studied  as  possible  starting 
points  for  the  application  of  the  dispersion  relationships  presented  in  this 
paper.  These  new  calculations  may  permit  the  use  of  a  longer  and  more 
realistic  mean  free  path.  Recently  other  models  have  been  developed  for  eQ 
that  could  help  moderate  the  height  and  broaden  the  resonance  with  larger 
values  for  the  mean  free  path.30 
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FIG.  1.  Real  (Nx)  and  imaginary  (N2)  parts  of  the  index  of  refraction  for 
Au/A1203(90%  A1203)  using  Drude's  model  with  an  effective  mean  free  path  of 
10A.  Dashed  line  represent  B  calculation,  solid-dashed  is  MG  and  dots  are 
experimental  data  from  H.  Craighead.21 

FIG.  2.  Changes  in  the  imaginary  part  of  the  dielectric  constant  due  to  the 
use  of  different  depolarization  factor  (D)  in  the  Maxwell  Garnett  equation  for 
the  dielectric  constant  of  Ag/MgO  (80%  MgO). 

FIG.  3.  Comparison  between  the  imaginary  part  of  the  dielectric  constant  as 
function  of  frequency  for  QSE  the  solid  line;  D  the  dashed  line;  and  KK  the 
sol  id- dashed  line. 

FIG.  4.  Comparison  between  the  experimentally  obtained  histogram  21  and  the 
log-normal  distribution  function  of  metallic  grains  radiuses  in  a  Ni/Al203 
(20%  A1203)  cermet  with  a  geometrical  standard  deviation  aln  of  1.3  and  an 
average  grain  size  of  30A. 

FIG.  5.  Unaveraged  and  averaged  QSE  dielectric  constant  of  Ag  metallic  grains 
for  a  grain  distribution  with  average  radius  of  30A  and  a  geometrical  standard 
deviation  of  1.3.  Solid  line  is  the  unaveraged  value  and  dashed  line  the 
averaged  one. 

FIG.  6.  Imaginary  part  of  the  index  of- refraction  for  Au/Al203  (90%  A1203). 
Interband  contribution  was  taken  from  H.  Ehrenreich  et  aL 18  Fig.  6(a) 
compares  BD,  BQ,  and  experimental  data.  21  Fig.  6(b)  compares  MGD,  MGQ,  and 
experimental  data.  Dotted  lines  are  used  for  D,  solid  line  for  QSE  modek  and 
dots  for  experimental  data. 

FIG.  7.  Comparison  between  experimental  data  21  and  calculations  for  N2.  (a) 
Ag/MgO  (80%Mg0)  cermet  compared  with  B  (dotted  line)  and  FSA  (solid  line) 
model  with  L  =  125A,  and  a  restricted  mean  free  path  of  10A.  (b)  Au/A1203 
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(84%A1*03)  cermet  compared  with  BO  (dashed  line),  BQ  (dotted  line)  and  FSAQ 
(solid  line)  models  with  L  =  80A,  a  =  30A,  and  restricted  mean  free  path  of 


FIG.  8.  Real  Nx  and  imaginary  part  N2  of  the  refractive  index  is  drawn  for  an 
Ag/MgO  (80%  MgO)  cermet.  Solid  line  is  the  RFSA  calculation,  dotted  lines 
represent  the  MG  model  and  dots  are  for  the  experimental  results.21 

(a)  Nt  and  N2  using  the  0  model  with  t  =  2.14  x  10  ies  in  both  MG  and  RFSA. 
The  radius  of  the  islands  is  taken  as  80A  and  L  =  6  islands. 

(b)  Nx  and  N2  using  QSE  with  a  restricted  mean  free  path  of  10A  for  dielectric 
constant  of  the  metallic  grains.  The  island  radius  are  taken  to  be  125A  and 
L  =  6  island. 

FIG.  9.  Real  and  imaginary  part  of  the  index  of  refraction  for  Ni/Al203, 
dotted  line  is  MGD,  solid  line  is  RFSA  and  dots  are  experimental  results.21 
(a)  61%  AL203  with  t  =  5.6  x  10  16s,  a  =  60A,  L  =  4  islands,  (b)  54%  A1203 
with  x  =  3.5  x  10  16s,  a  =  75A  and  L  =  4  islands. 
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